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COMPUTER  PROGRAM  FOR  CALCULATION  OF  SEPARATED 
' TURBULENT  FLOWS  ON  AXISYMMETRIC  AFTER- 
BODIES INCLUDING  EXHAUST  PLUME  EFFECTS 

by  Gary  D.  Kuhn. 

Nielsen  Engineering  & Research,  Inc. 
Mountain  View,  California 


1.0  INTRODUCTION 


One  of  the  critical  areas  in  the  design  of  both  aircraft 
and  missiles  is  the  interaction  between  propulsive  jets  and  the 
external  flow  over  the  aerodynamic  shapes  from  which  they  issue. 
The  drag  of  the  afterbody  and  exhaust  nozzle  can  be  a significant 
fraction  of  the  total  vehicle  drag.  In  order  to  obtain  know- 
ledge of  the  afterbody  and  nozzle  flow  fields  and  their  interac- 
tion as  early  as  possible  in  the  design  procedure,  it  is  desir- 
able to  have  a fast  computational  method  for  accounting  for  the 
viscous  flow  over  the  afterbody  including  separation  on  the  body 
and  .entrainment  of  the  flow  by  the  exhaust  jet.  Such  a method 
should  also  account  for  the  interaction  between  the  external 
inviscid  flow  and  the  viscous  flow. 

This  paper  describes  a method  for  predicting  the  viscous 
flow  field  about  an  axisymmetric  body  at  zero  angle  of  attack 
with  a cylindrical  sting  or  an  exhaust  jet.  Several  previous 
methods  for  calculating  such  flows  have  employed  empirical  or 
approximate  methods  of  predicting  the  location  of  separation 
and  the  characteristics  of  the  separated  flow  region  (c.f.  ref. 

1 and  2) . The  method  described  herein  uses  insights  from  those 
methods  to  establish  first  approximations  which  are  subsequently 
refined  by  iteration  to  satisfy  the  viscous-inviscid  interaction 
more  completely.  The  method  combines  a finite  difference  invis- 
cid-flow  method  with  integral  methods  for  the  boundary  layer  and 
exhaust  plume  mixing  layer.  An  eddy  viscosity  is  used  to  account 
for  the  Reynold's  stresses. 

The  basic  principle  of  the  viscous-inviscid  interaction 
calculation  is  to  define  a displacement  surface  for  the  viscous 
layer  which  forms  an  effective  boundary  for  the  inviscid  flow 
and  which  is  compatible  with  the  resulting  inviscid  pressure 
distribution.  When  the  boundary  layer  is  separated,  special 
techniques  are  needed  to  determine  the  appropriate  boundary 
layer  solution.  In  the  case  of  a body  with  a solid  sting,  bound- 
ary layer  theory  is  used  for  the  entire  viscous  region.  However, 
if  the  sting  is  replaced  by  an  exhaust  jet,  the  merging  of  the 
boundary  layer  with  the  exhaust  plume  mixing  layer  must  be 
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accounted  for.  Thus,  in  this  work,  an  entrainment  model  was 
developed  which  allows  the  viscous  interaction  between  the 
separated  boundary  layer  and  the  exhaust  plume  to  be  calculated. 

Since  the  principle  of  the  viscous-inviscid  interaction  is 
the  same  whether  a sting  or  a plume  exists  downstream  of  the  body, 
the  elements  of  the  calculative  method  will  be  described  in  the 
following  order.  First,  the  overall  viscous-inviscid  interaction 
method  will  be  described  in  order  to  provide  a framework  for 
explaining  the  various  individual  components  of  the  method. 

Next,  the  method  used  to  calculate  the  external  inviscid  flow 
will  be  presented.  That  presentation  will  be  brief  since  the 
method  has  been  described  in  detail  in  other  reports.  Following 
the  description  of  the  external  inviscid  flow  method  will  be 
the  method  used  to  calculate  the  inviscid  exhaust  plume  flow. 

Next,  the  boundary  layer  calculative  method  will  be  presented, 
followed  by  the  exhaust  plume  entrainment  model.  Finally,  some 
comparisons  with  experimental  data  are  shown,  illustrating  the 
capabilities  of  the  method.  The  remainder  of  this  report  con- 
tains a description  of  the  computer  programs  associated  with  the 
calculative  methods  followed  by  detailed  instructions  for  the 
use  of  the  programs.  The  composite  computer  program  described 
herein  contains  the  boundary- layer  and  entrainment  programs 
developed  for  this  work  and  an  inviscid-flow  program  which  is 
essentially  the  same  as  that  described  in  reference  3.  Both 
programs  are  incorporated  as  subprograms  into  a mainline  program 
which  performs  iterations  to  calculate  the  flow  for  axisymmetric 
bodies  of  the  type  used  for  nozzle  boattails. 

The  programs  are  written  in  the  FORTRAN  programming  language 
for  use  on  an  IBM  370  computer.  The  boundary-layer  program 
retains  the  capability  to  calculate  two-dimensional  flows  as  an 
option.  Both  the  boundary-layer  program  and  the  inviscid-flow 
program  can  be  used  individually,  without  iterating  if  so 
desired. 


2.0  VISCOUS-INVISCID  INTERACTION  METHOD 


Separation  is  the  result  of  an  adverse  pressure  gradient 
causing  reversal  of  the  low-energy,  low-momentum  fluid  near  the 
wall  in  the  boundary  layer.  The  adverse  pressure  gradient  is 
determined  by  the  inviscid  flow  field  and  is  the  net  result  of 
the  body  shape  and  the  displacement  effect  of  the  boundary  layer. 
Thus,  the  viscous  and  inviscid  flow  fields  are  interdependent. 
Accounting  for  this  interdependence  with  an  interaction  model  is 
made  difficult  by  the  fact  that  the  boundary-layer  equations  are 
parabolic  and  therefore  cannot  respond  to  disturbances  downstream 
of  local  stations,  while  for  subsonic  and  transonic  flows  the 
inviscid  flow  is  elliptical  and  therefore  subject  to  influence 
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by  the  entire  flow  field.  In  the  method  described  herein,  the 
viscous-flow  and  inviscid-flow  methods  are  used  alternately  in 
an  iterative  scheme. 

The  major  difficulty  encountered  in  calculating  viscous- 
inviscid  interactions  with  boundary  layer  separation  is  the  prob- 
lem of  obtaining  a smooth  surface  for  coupling  the  viscous  and 
inviscid  flows.  It  is  well  known  that  for  both  laminar  and  tur- 
bulent boundary  layers,  the  usual  boundary  layer  equations  exhib- 
it singular  behavior  at  the  point  of  zero  skin  friction  when  the 
pressure  gradient  is  prescribed.  The  details  of  the  method  used 
to  avoid  this  singular  behavior  will  be  described  in  Section  5 
of  this  report,  along  with  the  description  of  the  integral  bound- 
ary-layer method.  For  the  present  discussion,  it  is  sufficient 
to  note  that  a method  is  available  for  calculating  a smooth  sur- 
face which  represents  the  displacement  effect  of  the  boundary 
layer,  or  exhaust  plume  entrainment  layer,  including  boundary 
layer  separation.  The  surface  is  constructed  in  two  parts.  In 
regions  of  attached  flow,  or  far  downstream  in  the  exhaust  plume 
mixing  layer,  the  surface  shape  is  calculated  from  a solution  of 
the  viscous  flow  equations  with  the  velocity  at  the  edge  of  the 
boundary  layer  prescribed  from  a calculation  of  the  external 
inviscid  flow.  This  velocity  will  subsequently  be  called  the 
"inviscid  velocity."  In  regions  of  separated  flow,  the  edge 
velocity  is  calculated  from  the  viscous  flow  equations  with  the 
boundary  layer  displacement  thickness  prescribed.  This  velocity 
is  called  the  "viscous  velocity."  An  iteration  is  performed  to 
determine  the  surface  shape  so  that  the  "viscous  velocity"  matches 
the  "inviscid  velocity"  from  the  inviscid  flow  calculation. 

2.1  ESTIMATION  OF  BOUNDARY  OF  SEPARATED  REGION 

The  first  step  in  calculating  the  viscous-inviscid  interac- 
tion is  to  calculate  the  inviscid  flow  over  the  plain  body.  The 
resulting  distribution  of  the  velocity  at  the  boundary  is  then 
prescribed  as  the  boundary- layer- edge  velocity,  ue,  for  the 
boundary-layer  calculation.  If  strong  adverse  pressure  gradients 
exist,  the  boundary- layer  calculation  may  exhibit  a singularity 
where  the  skin-friction  coefficient  approaches  zero  and  the  nu- 
merical calculation  can  proceed  no  further  with  ue  prescribed. 
This  point  is  not  necessarily  the  true  location  of  separation, 
however,  since  the  interaction  with  the  inviscid  flow  has  not  yet 
been  accounted  for.  A first  approximation  to  the  separation 
point  location,  xs  is  obtained  by  a method  which  will  be  described 
in  Section  5,  along  with  the  detailed  description  of  the  boundary 
layer  method. 

At  the  separation  point,  xs  the  boundary-layer  calculation 
is  carried  on  into  the  separated  flow  using  a prescribed  distri- 
bution of  the  boundary  layer  displacement  thickness,  6*.  The 
result  of  that  calculation  is  a solution  for  the  boundary-layer- 
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edge  velocity  (the  "viscous  velocity")  which  may  or  may  not  agree 
with  the  "inviscid  velocity"  produced  by  the  inviscid  flow  theory. 
An  iterative  procedure  is  used  to  find  the  particular  variation 
of  6*  downstream  of  xs  for  which  the  "viscous  velocity"  and 
the  "inviscid  velocity"  agree. 

A two-parameter  analytical  formulation  is  used  to  represent 
the  effective  displacement  surface  between  separation  and  a point 
downstream  of  reattachment  or  in  the  plume  entrainment  region 
(fig.  1).  Thus,  the  effective  body  shape  is  given  by 


r = r,  + 5*  for 

w 

0 < x < x_ 
s 

and 

Xp  C X C 

CO 

(1) 

r ■ <rw  + S*>s  + (x  - 

xg) tan  6S 

for 

xs  < X < 

XP 

(2) 

A first  approximation 

to  the  angle , 

V 

is  obtained 

from 

an 

expression  presented  in  reference  4 . 

9 = tan  "*■  [ (r.,)  ] 

s w X 

+ 14.4  - 4.89Ms 

(degrees) 

(3) 

Downstream  of  Xp  the  boundary-layer-edge  velocity  is  again  pre- 
scribed on  the  viscous  flow  calculation  as  the  inviscid  boundary 
velocity  from  the  previous  calculation. 

2 . 2 ITERATION  PROCEDURE 

The  method  developed  in  this  work  consists  of  two  iterations. 
One  iteration  is  used  to  locate  the  separation  point,  xs  and 
determine  the  angle,  6S,  of  a conical  displacement  surface.  The 
second  iteration  calculates  the  best  solution  for  the  specific 
values  of  xs  and  0S,  The  iteration  procedure  is  described 
schematically  in  figure  2 . 

The  iteration  procedure  consists  of  two  cycles.  In  the  in- 
ner cycle,  the  inviscid  flow  and  the  boundary  layer  are  calculated 
alternately  until  the  largest  change  in  the  6*  solution  between 
iterations  becomes  smaller  than  a specified  percent  of  the  nozzle 
radius  and  simultaneously,  the  inviscid  calculation  procedure 
converges  to  a specified  tolerance.  At  each  step  of  the  cycle, 
the  boundary-layer  displacement  thickness  is  used  to  calculate 
an  augmented  body  shape  by  the  relation 

rn  ■ rw  + + (1  * (4) 

where  rn  is  the  effective  body  radius  at  iteration  n,  and  a 
is  a damping  factor,  usually  equal  to  0.5. 
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When  the  inner  cycle  has  been  terminated,  the  calculation  is 
complete  if  no  separation  occurred.  However,  if  separation  is 
present,  there  exist  two  solutions  for  the  boundary-layer-edge 
velocity  as  described  previously.  The  next  step  of  the  procedure 
is  to  calculate  the  squared  deviation  between  the  two  velocity 
solutions  downstream  of  at  corresponding  points.  The  accur- 

acy of  the  matching  of  the  two  solutions  is  then  indicated  by 
the  value  of  the  rms  error,  s,  where 


s 


N 

Z 

1 


ue 

v V 


u2  N 
eo  J 


1/2 


(5) 


where  ev  is  the  "viscous  velocity",  ej  is  the  "inviscid 
velocity",  and  ue0  is  the  free  stream  velocity.  This  quantity 
is  then  compared  with  the  value  from  the  previous  iteration  and 
if  a minimum  has  not  been  found,  xs  and  0S  are  adjusted  in  an 
appropriate  manner,  the  boundary  layer  is  recalculated  using  the 
new  values  of  x?  and  0S  and  the  calculation  reenters  the 
inner  cycle.  This  process  is  repeated  until  a minimum  of  the 
rms  error,  s,  is  found  or  until  the  value  of  s is  less  than  a 
specified  amount.  For  practical  purposes,  a value  of  s of  0.01 
or  less  is  usually  adequate. 


The  manner  in  which  xs  and  0S  are  adjusted  will  now  be 
described.  The  first  value  of  s (equation  (5))  is  calculated 
using  the  estimated  values  of  x-  and  9g.  The  next  step  is  to 
change  0S  a small  amount  (.25  degree)  and  calculate  a second 
value  of  s.  The  third  step  is  to  change  a small  amount, 

keeping  0S  at  the  original  value,  to  calculate  a third  value 
of  s.  With  these  three  solutions,  the  gradient  of  s with 
respect  to  xg  and  0S  is  calculated  from 


0s  s2~sl 
S0  - 30s  " A0s 


(6) 


s 


X 


(7) 


The  maximum  rate  of  decrease  of  s with  respect  to  and 

0g  should  be  in  the  direction  of  the  negative  of  the  'gradient. 
Therefore,  both  xs  and  0S  are  adjusted  according  to 


Ax  = -k  (sin  a | L 

s sxl 


(8) 
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Sg 

A 0 = -k  cos  a 

s ls0 


(9) 


where 

k = .003 

L = the  length  of  the  body 

o-*-’1  (=?*■] 

After  the  third  iteration,  both  quantities  are  adjusted  along  the 
direction  of  the  negative  gradient  until  either  the  tolerance  is 
met  or  a minimum  value  of  s is  reached.  If  a minimum  is  found 
without  satisfying  the  tolerance,  the  sequence  is  restarted, 
calculating  a new  gradient  direction,  using  half  steps  for  each 
subsequent  adjustment  of  x and  0S.  If  this  procedure  does 
not  produce  a value  of  s less  than  the  specified  tolerance 
after  four  restarts,  the  sequence  is  terminated. 

A third  element  in  the  iteration  procedure  is  the  calcula- 
tion of  the  inviscid  exhaust  plume.  In  the  theory  to  be  presen- 
ted subsequently,  the  entrainment  of  the  boundary  layer  flow  and 
the  turbulent  mixing  between  the  external  inviscid  flow  and  the 
exhaust  plume  flow  are  assumed  to  take  place  along  a boundary 
defined  by  the  radius  of  the  exhaust  plume  with  all  mixing  ne- 
glected. This  "inviscid  reference  surface"  is  calculated  twice 
during  the  iteration  process  described  previously.  The  first 
plume  calculation  is  performed  assuming  a constant  external 
pressure.  This  gives  a first  approximate  jet  boundary  upon  which 
the  viscous-inviscid  iteration  can  be  initiated.  After  4 cycles 
of  the  iteration,  a new  exhaust  plume  is  calculated  using  the 
external  pressure  distribution  that  exists  at  that  time  as  cal- 
culated by  the  inviscid  theory.  Tests  in  which  the  plume  was 
calculated  a third  time  using  the  final  converged  pressure  dis- 
tribution indicated  that  only  a slight  change  in  the  plume 
occurred.  The  effect  on  the  flow  characteristics  on  the  body  was 
small.  Thus,  the  solution  is  not  sensitive  to  the  exact  shape  of 
the  jet  mixing  boundary. 


3.0  INVISCID-FLOW  MODEL 


The  inviscid-f low  method  employed  herein  is  the  method  of 
South  and  Jameson  described  in  reference  5 . The  computer  program 
used  is  that  described  in  reference  3 with  minor  modifications 
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/■ 


to  accommodate  the  iteration  between  the  inviscid  and  viscous 
flows.  The  detailed  derivation  of  the  flow  equations,  the  cal- 
culative  method,  and  the  computer  program  are  contained  in  ref- 
erences 3 and  5 and  will  not  be  repeated  here. 

The  inviscid-f low  method  calculates  the  flow  field  about 
axisymmetric  bodies  by  a finite-difference  solution  of  the  exact 
equation  for  the  disturbance  potential  for  axisymmetric  compress- 
ible flow.  The  only  simplifying  assumption  that  is  used  is  that 
the  flow  is  isentropic.  Thus,  the  entire  subsonic  velocity  range 
can  be  calculated,  and  slightly  supersonic  free-stream  Mach 
numbers  are  permissible  so  long  as  entropy  production  by  shock 
waves  is  negligible.  In  addition,  blunt  as  well  as  slender  bodies 
can  be  treated  accurately.  A general  description  of  the  program 
is  summarized  here  from  reference  3. 

"One  of  the  important  considerations  when  trying  to  solve 
the  full  potential  equation  is  the  choice  of  a coordinate  system. 
For  complex  three-dimensional  shapes  cartesian  coordinates  may 
be  best;  however,  for  simpler  two-dimensional  or  axisymmetric 
shapes  the  use  of  a coordinate  transformation  such  that  the 
body  lies  along  a coordinate  line  can  greatly  simplify  the  appli- 
cation of  the  exact  boundary  condition  at  the  body  surface.  The 
program  described  in  this  paper  uses  a body-normal  coordinate 
system  for  closed  bodies.  For  open  bodies  (i.e.,  bodies  with  a 
sting  or  simulated  wake)  it  uses  a body-normal  system  on  the 
forebody  up  to  the  first  horizontal  tangent  and  a sheared  cylin- 
drical coordinate  system*  aft  of  that  point.  This  coordinate 
system  is  suitable  for  closed  bodies  which  are  blunt  on  both  ends 
and  convex  and  smooth  over  the  entire  body  or  for  open  bodies 
which  are  blunt-nosed  and  convex  and  smooth  up  to  the  first  hori- 
zontal tangent.  It  is  possible  to  treat  pointed  bodies  and 
bodies  with  slope  discontinuities  but  the  coordinate  system  is 
not  well  suited  for  them  and  their  solution  may  not  be  so  accur- 
ate as  the  blunt-body  solutions . " 

"A  stretching  is  applied  to  both  the  normal  and  tangential 
coordinates  such  that  the  infinite  physical  space  is  mapped  to 
a finite  computational  space.  Thus,  the  boundary  condition  at 
infinity  can  be  applied  directly  and  there  is  no  need  for  an 
asymptotic  far-field  solution.  Details  about  the  stretching 
functions  are  given  in  (reference  3)  appendix  A." 

"The  general  method  of  solution  is  to  replace  the  governing 
second-order  partial  differential  equation  with  a system  of 


The  origin  of  the  x,  r,  0 coordinate  system  is  shifted  from 
the  body  axis  to  the  body  surface  in  a parallel  transformation. 
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finite-difference  equations,  including  Jameson's  rotated  differ- 
ence scheme  at  supersonic  points.  The  difference  equations  are 
solved  by  a column  relaxation  method." 

"The  boundary  condition  at  the  body  surface  is  applied 
through  the  use  of  dummy  points  inside  the  body.  Details  of 
this  computation  are  given  in  (reference  3)  appendix  B."  / 


4.0  INVISCID  PLUME  MODEL 


The  inviscid  exhaust  flow  calculation  method  used  in  this 
work  is  the  same  as  the  Henson  and  Robertson  shock -expans ion/one- 
dimensional method  (reference  6)  with  two  exceptions: 

1.  Compressions  are  computed  as  reverse  Prandtl-Meyer  ex- 
pansions rather  than  by  means  of  Newtonian  impact  theory. 

2.  Flow  areas  are  computed  as  plane  areas  normal  to  the 
nozzle  centerline  rather  than  as  spherical  areas  (normal 
to  the  nozzle  centerline  and  the  local  plume  contour) . 

A diagram  of  the  nozzle  plume  region  is  shown  in  figure  3.  The 
plume  shape  is  calculated  in  a series  of  straight  line  segments 
which  approximate  the  curved  shape  of  the  plume.  The  external 
pressure  is  compared  with  the  internal  pressure  at  the  exit  plume 
boundary  point  and  at  each  consecutive  boundary  point,  and  an  ex- 
pansion or  compression  turn  is  calculated  depending  on  the  com- 
parison. 

The  expansion  or  compression  at  the  nozzle  exit  is  computed 
by  means  of  the  Prandtl-Meyer  relationships  from  specified  values 
of  the  nozzle  expansion  angle  and  the  exhaust  jet  Mach  number  such 
that  the  local  external  pressure  is  matched.  This  calculation 
yields  M.-,  t^ie  Mac^  number  at  the  upstream  end  of  the  first 
plume  segment  and  the  inclination  of  this  segment  (the  local  flow 
direction) . The  local  flow  angle  is  used  to  compute  the  radius 
at  the  next  downstream  point,  r2  from  which  follows  the  area,  A2 . 
The  Mach  number  at  the  downstream  end  of  the  segment,  M2^,  is 
then  computed  from  the  one-dimensional  relationship 

(A/A* ) 21  = (A/A*)12  a2/Ax  (10) 

The  Mach  number  at  the  upstream  end  of  the  second  segment,  M-2, 
is  then  determined  by  matching  the  internal  static  pressure  to 
the  local  external  static  pressure  at  point  2 by  means  of  an 
expansion  or  compression  beginning  at  a flow  direction  equal  to 
that  of  the  first  segment  and  the  Mach  number  M2j.  The  Mach 
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number  at  the  downstream  end  of  the  segment,  , is  then  computed 
by  applying  equation  (10)  to  the  second  segment.  Successive  seg- 
ments are  constructed  in  the  same  manner  until  the  plume  boundary 
becomes  parallel  with  the  nozzle  centerline. 


The  plume  boundary  velocity,  q , corresponding  to  each  axial 
station  is  computed  from  isentropic*  relations  with  the  specified 
nozzle  stagnation  conditions.  Thus 

q = a.M0- 

t 22 

v 


1 + §<Yp  - DM2 


22 


-1/2 


(11) 


where 


at  = The  sound  speed  corresponding  to  the  nozzle 
stagnation  temperature 


= The  ratio  of  specific  heats  of  the  plume  gas 


In  order  to  approximate  the  velocity  distribution  in  each  cross 
section  of  the  plume,  the  flow  is  assumed  to  be  locally  source- 
like, so  that  the  velocity  components  required  at  r , the  inner 
boundary  of  the  turbulent  mixing  layer  are  approximated  by 


u 


c 


v 

c 


r r 

dr  i 

2 

1 4- 

c 

w 

■L  r 

r 

dx 

^ w 

/ 

-1/2 


u 


c r. 


w 


ffw 

dx 


(12) 

(13) 


5.0  DEVELOPMENT  OF  BOUNDARY-LAYER  METHOD 


The  complete  derivation  of  the  governing  equations  for  the 
boundary  layer  has  been  presented  in  reference  7.  In  this  report 
the  derivation  is  summarized  and  modified  for  the  application  of 
interest  herein. 

5.1  ASSUMPTIONS 

The  analysis  is  based  on  the  following  assumptions: 

1.  The  governing  equations  are  those  for  a compressible 
turbulent  boundary  layer. 

2.  The  air  behaves  as  an  ideal  gas. 

3.  The  molecular  viscosity,  y,  is  proportional  to  the 
temperature. 
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4.  The  specific  heat  of  the  gas  is  constant. 

5.  The  wall  is  either  two  dimensional  or  axisymmetric,  but 
can  have  an  arbitrary  profile  in  the  direction  of  flow 
as  long  as  the  longitudinal  radius  of  curvature  of  the 
wall  is  large  compared  to  the  boundary  layer. 

6.  The  pressure  is  constant  through  the  boundary  layer 
normal  to  the  wall. 

7.  The  wall  temperature  is  constant. 

5.2  BOUNDARY- LAYER  EQUATIONS  FOR 

COMPRESSIBLE  TURBULENT  FLOW 

The  basic  notation  and  coordinate  scheme  are  shown  in  fig- 
ure 4.  Note  that  the  same  symbols  are  used  for  the  physical 
coordinates  of  both  two-dimensional  and  axisymmetric  configura- 
tions. Thus,  r denotes  the  distance  of  a point  from  the  axis 
of  an  axisymmetric  configuration,  x is  the  distance  along  the 
axis,  measured  from  the  nose  and  the  dimension  y is  measured 
from  the  body  surface  normal  to  the  axis.  It  is  noted  that  these 
coordinates  are  not  the  usual  boundary  layer  coordinates . They 
are  employed  to  avoid  difficulties  at  sharp  concave  corners  such 
as  a boattail-sting  junction. 

The  governing  equations  describing  the  steady  flow  of  a 
compressible  turbulent  boundary  layer  in  these  coordinates  are: 

(rpu)  + {rp [v  -(r  ) u] } = 0 (14) 

a w x y 

puux  + p[v  - (rw)  ]Uy  = -px  + (1/r)  (rue  uy)y  (15) 

puS  + p [v  -(r  ) u]  S = (1/r)  (ry£S  ) (16) 

x w x y y y 

where  p is  the  density,  u and  v are  the  x and  y velocity 
components,  3 is  the  eddy  viscosity  factor,  y is  the  molecular 
viscosity,  p is  the  pressure  and 

S = Tt/Tt  - 1 (17) 

e 

where  Tfc  is  the  total  temperature  and  e denotes  the  boundary 
layer  edge.  Both  laminar  and  turbulent  Prandtl  numbers  are 
assumed  to  be  unity. 

Equations  (14)  , (15)  , and  (16)  are  easily  applicable  to 
laminar,  transitional  and  turbulent  flow.  In  laminar  flow, 
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substitution  of  B = 1 reduces  the  equations  to  those  for  a 
laminar  boundary  layer.  Further,  suitable  variation  of  the  eddy 
viscosity  makes  the  equations  applicable  to  the  transition  region. 

5.3  TRANSFORMATION  OF  AXISYMMETRIC 
BOUNDARY-LAYER  EQUATIONS 

The  Probstein-Elliott  transformation  is  (ref.  8) 

dx  = [rw(x)/L]2dx  (18) 

dy  = [r(x,y)/L]dy  (19) 

where  L is  an  arbitrary  reference  length,  rw(x)  is  specified 
by  the  body  shape  and  r(x,y)  is  given  by 

r (x,y ) = rw(x)  + y (20) 

The  transformed  continuity  equation  has  the  form 

(pu)~  + (pv)~  = 0 (21)  ^ 

where 

pu  = pu  (22) 

v 

and 

* \ f -s 

pv  = rL/r2  p v- (r  ) u + L2/r2  y pu  (23) 

. WJ  x 

Applying  the  transformation  to  the  momentum  and  energy  equations 
(15)  and  (16)  yields  the  transformed  equations 

puu~  + pvu~  = -p~  + [ (1 + kty)uBu~l - (24) 

puS~  + pvS~  = [ (1 + kty)yBS~]~  (25) 

a y y y 

where  t is  the  transverse  curvature  factor 

t = 2L/r2  (26) 

and,  from  equations  (19)  and  (20) , 

y = (rw/L)y  + y2/2L  (27) 
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For  flows  in  which  the  transverse  curvature  is  small,  letting 
k = 0 in  equations  (24)  and  (25)  produces  the  equations  of  a, 
two-dimensional  boundary  layter . The  transverse  curvature  terms 
may  be  negligible  for  an  axisymmetric  flow  if  the  body  radius  is 
large  compared  to  the  boundary- layer  thickness.  The  Probstein- 
Elliott  transformation  is  thus  a first-order  correction  of  the 
approximate  equations  for  the  effect  of  transverse  curvature, 
allowing  the  boundary-layer  thickness  to  be  of  the  same  order  as 
the  body  radius. 

5.4  TRANSFORMATION  OF  THE  COMPRESSIBLE 
BOUNDARY-LAYER  EQUATIONS 


The  Stewartson  transformation  (ref.  9)  along  with  the 
assumption  that  the  viscosity  varies  linearly  with  the  tempera- 
ture and  laminar  and  turbulent  Prandtl  numbers  are  unity,  re- 
duces the  equations  of  the  compressible  boundary  layer  to  those 
of  an  incompressible  boundary  layer. 

Ux  + VY  = 0 (28) 

UUY  + VUv  = (S  + 1)U  U + CV  [ (1  + kty) 0U  ] (29) 

X Y e eX  eo  Y Y 

USX  + VSy  = Cve  [ (1  + kty) esy]  (30) 

O Y 

where  the  subscript  o refers  to  reference  conditions  in  the 
undisturbed  free  stream. 


The  coordinate  y is  not  transformed  in  the  transverse 
curvature  terms  because  integration  of  the  equations  across  the 
boundary  layer  is  anticipated  and  only  corresponding  values  are 
needed  in  those  terms.  The  Chapman-Rubesin  parameter,  C is  a 
constant  evaluated  at  the  wall  temperature;  that  is. 


C = 


u T 
pw  e 


/ 

° i l 


H = (v*.„r  (T< 


+ T. 


■/ 

fT  + T 

w s 

. . 

(31) 


where  Sutherland's  law  is  used  to  evaluate  the  viscosity,  with 
Ts  a constant. 

In  the  remainder  of  this  report,  the  solution  of  the  energy 
equation  (30)  will  be  approximated  for  boundary  layers  by  the 
Crocco  relation 

S = Sw(l  - U/Ue)  (32) 
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The  density  profiles  are  then  related  to  the  velocity  profiles 
through  the  temperature. 

Pe/P  = T/Te  = Tt/Te  - u2/(2gJcpTe)  (33) 

but 

cp  = yR/(y-l)  (34) 

So 

T/Te  = 


or 

(Pe/P)  = (T/Te)  = (S  + 1) (1  + m£)  - me(u/ue)2  (36) 

where 

me  = 2 (37) 

Thus,  the  velocity  profiles  found  to  be  valid  for  incom- 
pressible, two-dimensional  turbulent  boundary  layers  can  be  used 
by  simply  transforming  the  input  quantities  to  the  incompressible 
plane,  performing  the  calculation  for  an  equivalent  incompressible 
boundary  layer,  and  then  transforming  the  results  back  to  the 
compressible  plane,  and  for  an  axisymmetric  flow,  back  to  the 
axisymmetric  coordinates. 

5.5  DEVELOPMENT  OF  INTEGRAL  BOUNDARY -LAYER  METHOD 
5.5.1  Integral  Equations 

The  integral  method  used  herein  was  described  in  detail  in 
reference  10.  Families  of  integral  equations  are  derived  by 
eliminating  V between  the  momentum  and  continuity  equations  (28) 
and  (29)  , and  then  taking  weighted  integrals  of  the  resulting 
equation  across  the  boundary  layer. 

r&  . 

J - UY 

o o 


uxdn  - 


(S  + 


l>Wx 


VTt. 


Ft  /Tj 


- uV(2gJYRTe) 


(35) 


- Cv  [(1  + kty) 3uv]v}f (Y)dY  = 0 (38) 

o 
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In  the  present  case,  the  functions 

f = Yn  ; n = 0,  1 (39) 

produce  the  momentum  and  moment  of  momentum  integral  equations, 
respectively. 

5.5.2  Velocity  Profiles 

The  Y dependence  of  the  integral  equations  is  eliminated 
by  substituting  an  appropriate  parametric  formulation  for  the 
velocity  profiles.  The  function  used  for  the  present  theory  is 
a modification  of  the  Coles  family  of  profiles  (ref.  11)  with 
a laminar  sublayer  added  and  the  wake  function  approximated 
analytically . ■ 

U - Ux[2.5  In (1  + Y+)  + 5.1  - (3.39Y+  + 5 . 1 ) exp (-0 . 37Y+) ] 

+ V1"2  7 T[  (40) 

The  parameter  UT  is  the  friction  velocity 

UT  = <Cf/|cf  [)Ue/|Cf  |/2  (41) 

The  variable  Y+  is  defined  to  account  for  the  axisymmetry  of 
the  flow 

Y+  = (L/rw) (|UT|Y/v  ) (42) 

o 

The  other  parameters  in  equation  (40)  are  6^,  the  transformed 
boundary-layer  thickness  and  Ug  a wake  velocity.  The  expon- 
ential terms  and  the  additional  unit  in  the  logarithmic  term 
provide  a smooth  transition  from  the  turbulent  flow  to  the  wall 
through  a laminar  sublayer. 

5.5.3  Eddy  Viscosity 

The  eddy-viscosity  model  used  in  this  work  is  an  extension 
of  the  two-layer  model  used  by  Kuhn  (ref.  7)  including  an  inter- 
mittency  function  for  the  outer  layer  and  a modification  of  the 
outer  layer  for  adverse  pressure  gradients  and  separated  flows. 

In  the  inner  layer  of  attached  flows,  the  eddy-viscosity  param- 
eter, 3,  is  represented  by  an  exponential  expression  based  on 
the  law  of  the  wall.  In  the  outer  layer  Clauser's  expression, 
modified  for  adverse  pressure  gradients,  is  used  along  with  an 
intermittency  function  giving 

3 = [0.013+0.0038  exp(-6*px/15Tw) ]Ue6j/[l+5.5(y/6)6}  (43) 
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where 

* 

«Sk  = I (1  - u/ue)dy  (44) 

o 

For  favorable  pressure  gradients,  the  exponential  term  in  equa- 
tion (43)  is  taken  to  be  unity. 

For  separated  flows,  the  eddy  viscosity  across  the  entire 
layer  is  represented  by  a relation  based  on  the  velocity  profile 
above  the  U = 0 line. 

6 ~1  f ^ 

8 = 0.01311  + 5.5(y/«rj  (Ue/v)l  (1  - u/ue)dy  (45) 

^u=0 


5.5.4  Transitional  Eddy  Viscosity 

Transition  from  laminar  to  turbulent  flow  is  calculated  by 
letting  the  eddy  viscosity  change  from  a laminar  viscosity  to  a 
fully  turbulent  value  over  a short  distance  according  to  the 
relation  (ref.  7) 

Bt  * (1  - exp [-K (x  - xt)2])(6T  - 1)  + 1 (46) 

where  < has  been  chosen  to  provide  a transition  length  of 
approximately  ten  boundary  layer  thicknesses,  so  that 

k = 0.0001/62 

and  is  the  value  of  6 at  the  beginning  of  transition, 

Xt*  The  location  of  Xfc  must  be  specified  independently. 

5.5.5  Equations  Solved  f 

Substitution  of  equation  (40)  into  the  two  equations  pro- 
duced by  equations  (38)  and  (39)  produces  two  ordinary  differ- 

UT'  V 
evaluating 
Ug  from 

(47) 


(48) 


ential  equations  for  the  variation  of  the  variables 
<5  . , and  U with  x.  A third  equation  produced  by 
equation  (50)  at  Y = 6 allows  the  elimination  of 
the  equations,  leaving  a set  of  two  equations 


\l(Vx  + \2Si  + A13<Vx“Bn 


A21(Vx  + A225i  * ^jl'e'x  = '! 
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where  the  coefficients  are  integrals  evaluated  from  the  follow 
ing  definitions. 


For  n = 1 or  2 


6. 


Gn  ■ 2UY11-1  - (n  - l)j  UdY  - fij"1  + (n  - l)!1  UdY 


o 


I = 
n 


$i  n 

1 Yn  ± (S  + l)dY 


= -CUt|ut|k 


(49) 

(50) 

(51) 


(1  + kty)e||dY 


(52) 


with 


K = 


w 


1 + m 


1 - 3y 
' 2 (Y-l) 


1 + m 


o' 


The  Anj  and  Bn  of  equations  (47)  and  (48)  become 


or  2 , 

_ 'ii 

H 

i 

o 

?i 

•"N 


G dY 

nan 

T 


An2  “ j GnW7  dY 

i l 


o 
6 . 


A - = f G dY  - U I 

n3  I n3U  e n 

* P 


B = J - (r 
n n 


W>x  j1  Gn^t  dY 


J 


(53) 


(54) 
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These  integrals  are  evaluated  numerically,  using  a simple  11- 
point  trapezoidal  integration.  More  sophisticated  approaches 
were  examined,  but  were  not  found  to  effect  the  calculations 
significantly.  The  simpler  method  was  therefore  chosen  in  the 
interest  of  speed  and  cost. 


The  coefficients  A^  and  the  term  Bn  are  functions  of 
the  variables  UT , 62,  ana  Ue.  The  usual  procedure  for  solving 
equations  (47)  and  (48)  for  attached  boundary  layers  is  to  pre- 
scribe the  pressure  distribution  or  the  boundary- layer -edge 
velocity  distribution,  UQ.  However,  if  separation  occurs,  the 
pressure  distribution  cannot  be  prescribed  arbitrarily  in  the 
separated  region.  If  an  adverse  pressure  gradient  is  prescribed 
for  an  attached  boundary  layer,  the  value  of  UT  can  approach 
zero.  When  UT  vanishes,  the  coefficients  A]^  and  A22  in 
equations  (47)  and  (48)  also  vanish,  producing  a singularity. 

The  singularity  can  be  removed  by  rearranging  the  equations  so 
that  UT  is  not  a dependent  variable.  One  method  of  accomplish- 
ing this  is  simply  to  rearrange  the  equations  so  that  UT  can 
be  prescribed  and  Ue  can  be  calculated  as  a dependent  variable 
as  shown  in  reference  10. 

A different  method  of  avoiding  the  singularity  at  UT  = 0 
is  used  herein.  The  displacement  thickness;  6*,  is  expressed  in 
terms  of  UT , 5^,  and  Ue  through  the  definition 


6*  = J (1  - pu/peue) (r/rw) dy  (55) 

o 


The  result  is  differentiated  with  respect  to  x,  producing  a 
third  equation. 


A3l'Vx  + A326i  + A33(Vx  + A34^  = B3 

3C 


(56) 


where  the  A^^  and  E-  are  functions  of  UT , 6^,  and 
evaluated  from  the  following  definitions 


U 


O' 


X3  " 


U, 


G - _L  ffel  _ _1_ 
G3  3U  [ P J Ue 


f* 

a 

N + il 

J 

1 Ue 

p J D 2> 

_1_ 

u„ 


Y + i Me 


Kfi’ 


(57) 

(58) 
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and 


<5*  . , 

— (r»> 
w x 


1 + m 


1 + m 


eo> 


y + 1 

2 (y-1) 


(59) 


With  these  definitions,  the  coefficients  A3j 
are  evaluated  by  equations  (54)  for  j = 1,2, 
coefficient,  A34  is  given  by 

Y + 1 


A34  = 


_w 

L 


1 + m 


1 + m. 


2 (y-1) 


and  the  term  B3 
and  3.  The  final 


(60) 


Like  the  coefficients  A33  and  A21,  the  coefficient  A31 
also  vanishes  when  U =0.  However,  it  can  be  shown  that 
A]^/A3i  and  ^2i/A3iT  are  both  finite  when  UT  = 0.  This  allows 
the  three  equations,  (47),  (48),  and  (56),  to  be  reduced  to  two 
ordinary  differential  equations  in  the  four  dependent  variables , 

UT , 6*,  5^,  and  Ue.  The  resulting  equations  are: 

all6ix  + a12<ue>  = bl  <61> 

3v 


where,  with  n 


and 


a216ix  + a22 (Ue3  x “ b2 
1 and  2 and  j = 1 and  2 

anj  = (Anl/A31)A3, j+1  “ An,j+1 


bn  * (Anl/A31} (B3 


- U 6* ) -B 
ex  n 


(62) 


(63) 


(64) 


The  value  of  UT  can  be  obtained  directly  from  known  values  of 
6j_,  5*  and  Ue  by  solving  the  nonlinear  relation  (55) . 

Prescribing  the  distribution  of  6*  in  equations  (61)  and 
(62)  and  obtaining  UT  from  equation  (55)  is  equivalent  to  pre- 
scribing UT  as  described  previously. 

5.5.6  Method  of  Solution  of  Equations 


The  equations  (47)  and  (48)  , or  (61)  arid  (62)  are  integrated 
numerically  using  the  following  procedure:  If  <|>  denotes  any 
one  of  the  variables  6^,  Ue,  or  UT  , the  value  of  dtf>/dx  at  a 
particular  x-station  can  be  calculated  from  the  equations  as 
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functions  of  the  variables  themselves.  The  value  of  <(>  at  the 
next  x- station  is  then  found  by  means  of  the  elementary  predic- 
tor-corrector scheme  described  by  Nash  and  Hicks  (ref.  12) . Thus, 

’"x+ix/a  = + <a<e/dx>x  • ¥ 1 

> (65) 

♦x+Ax  " *x  + <d^dx)x+Ax/2  ' ix  J 

In  the  present  work,  Ax  is  taken  to  be  l/16th  of  each  finite- 
difference  step  of  the  inviscid  calculational  mesh.  This  pro- 
vides a fast,  accurate  method  of  solution. 

5.6  CALCULATION  OF  SEPARATED  BOUNDARY  LAYERS 

As  discussed  previously,  it  is  well  known  that  for  both 
laminar  and  turbulent  boundary  layers , the  usual  boundary  layer 
equations  exhibit  singular  behavior  at  the  point  of  zero  skin 
friction  when  the  pressure  is  prescribed  (or,  equivalently,  when 
the  boundary- layer- edge  velocity  is  prescribed) . Stewartson 
(ref.  13)  suggests  that  if  the  boundary  layer  is  to  continue 
downstream  of  separation,  the  external  stream  must  adjust  itself 
so  that  the  singularity  cannot  appear.  Furthermore,  he  says 
that  the  solution  downstream  of  separation  is  no  longer  specified 
uniquely  by  the  boundary-layer-edge  velocity  and  the  velocity 
profile  at  the  beginning  of  the  boundary  layer.  This  is  another 
way  of  saying  that  the  boundary-layer-edge  velocity  cannot  be 
specified  arbitrarily  at  the  separation  point,  or  in  the  separa- 
ted flow  region  downstream.  If  the  solution  continued  downstream 
of  separation  does  not  satisfy  some  downstream  boundary  condition, 
the  separation  point  must  move  until  the  downstream  condition  is 
met.  It  has  been  observed  by  many  authors  (c.f.  Cebeci  et  al, 
ref.  14  and  Gerhart  and  Bober,  ref.  15)  that  the  use  of  the  ex- 
perimental pressure  distribution  from  a separated  flow  in  a solu- 
tion of  the  boundary  layer  equations  does  not  always  produce  a 
singularity.  Often  the  skin  friction  will  simply  decrease  to 
some  minimum  value  and  then  increase,  with  the  calculated  boundary 
layer  remaining  attached  in  a region  known  from  experiment  to  be 
separated!  The  reason  for  this  is  believed  to  be  due  to  the 
neglect  of  the  terms  in  the  Navier-Stokes  equations  which  respond 
to  downstream  perturbations,  i.e.,  the  normal  stress  terms  and 
normal  pressure  gradient  terms.  On  the  other  hand,  the  work  of 
Newman  (ref.  16)  and  Simpson  (ref.  17)  has  shown  that  these  terms 
are  negligible  a short  distance  away  from  separation  even  though 
they  can  be  very  significant  in  the  immediate  vicinity  of  the 
separation  point. 

The  approach  taken  in  the  present  work  is  an  engineering 
approximation  based  on  the  following  assumptions: 
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1.  Upstream  influences  are  transmitted  predominantly  through 
the  external  inviscid  flow  which  is  elliptical  in  nature 
for  a locally  subsonic  flow. 

2.  Violations  of  the  boundary  layer  assumptions  such  as 
significant  normal  pressure  gradients  associated  with 
separation  can  be  neglected  on  the  grounds  that  they  are 
only  important  in  a small  neighborhood  of  the  separation 
point  and  have  negligible  effect  on  the  rest  of  the  flow. 

3.  The  effective  boundary  which  determines  the  external  flow 
field  is  a smooth  boundary.  That  is,  discontinuities 

in  body  slope  are  smoothed  out  by  the  boundary  layer. 

The  significance  of  assumption  number  3 is  that  a calculative 
technique  which  allows  singular  behavior  at  the  separation  point 
cannot  produce  a smooth  body  shape  and  may  lead  to  instabilities 
in  the  iteration  between  the  viscous  and  inviscid  flow  calcula- 
tions. The  approach  taken  herein  is  to  obtain  information  about 
the  location  of  the  separation  point  independently  from  the 
boundary  layer  calculation,  and  simply  let  the  skin  friction  be 
zero  at  that  point  keeping  the  boundary- layer-edge  velocity  and 
the  displacement  thickness,  6*  continuous.  This  means  that  the 
variation  of  upstream  of  the  separation  point  is  probably 

slightly  in  error,  the  amount  depending  upon  the  extent  of  the 
neglected  upstream  influence.  Also,  other  important  boundary 
layer  quantities  may  be  slightly  in  error  in  that  region.  For 
many  practical  cases,  however,  the  errors  are  not  expected  to 
be  significant. 

As  discussed  previously,  if  strong  adverse  pressure  gradients 
are  prescribed  in  equations  (47)  and  (48) , the  boundary- layer 
calculation  may  reach  a point  where  the  friction  velocity,  UT 
approaches  zero  and  the  numerical  calculation  can  proceed  no 
further  with  ue  prescribed.  This  point  is  not  necessarily  the 
true  location'  of  separation,  since  the  interaction  with  the  in- 
viscid flow  has  not  been  accounted  for. 

A better  approximation  to  the  location  of  separation  xg 
can  be  found  by  computing  the  shape  factor,  Hi . 

6i  6i 

Hi  = | (l-U/Ue)dy/j  U/Ue(l-U/Ue)dy  (66) 

o o 

For  the  velocity  profiles  given  by  equation  (40)  this  parameter 
has  a value  of  4.0  when  UT  = 0.  However,  experimental  measure- 
ments indicate  that  separation  actually  occurs  when  Hi  is 
approximately  2.0.  This  suggests  that  the  present  boundary -layer 
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formulation  is  not  accurate  in  the  vicinity  of  the  separation 
point.  This  is  not  unexpected  in  the  light  of  the  previous 
discussion.  The  velocity  profiles  of  equation  (40)  have  been 
demonstrated  to  provide  a reasonable  approximation  to  both  at- 
tached and  separated  boundary  layers  (ref.  18)  but  to  be  of 
questionable  accuracy  at  the  point  of  zero  skin  friction.  How- 
ever, the  effect  of  this  inaccuracy  is  believed  to  be  confined 
to  the  immediate  vicinity  of  the  separation  point  for  the  kind 
of  flows  of  interest  here.  In  order  to  obtain  a smooth  solution 
the  first  approximation  to  the  value  of  xs  is  taken  to  be  the 
location  at  which  Hi  equals  a specified  value  for  the  inviscid 
velocity  calculated  for  the  plain  body. 

When  a value  of  xs  is  established,  the  boundary  layer  is 
calculated  frpm  the  nose  of  the  body  to  xs  with  the  distribu- 
tion of  ue  'from  the  external  inviscid  flow  calculation  pre- 
scribed using  equations  (47)  and  (48) . At  xs  the  velocity  ue 
and  the  displacement  thickness,  6*  are  assumed  to  be  continuous 
(see  figure  5).  The  friction  velocity,  UT , on  the  other  hand,  is 
assumed  to  be  zero  as  the  initial  value  for  the  integration  into 
the  separated  region  using  equations  (61)  and  (62)  with  6* 
prescribed.  It  must  be  noted  that  the  calculation  could  be 
started  with  a positive  value  of  UT  and  the  calculation  would 
proceed  through  the  point  where  UT  = 0 for  an  appropriate  6* 
distribution  with  no  singularity.  However,  due  to  the  approxi- 
mations inherent  in  the  theory  in  the  vicinity  of  the  separation 
point,  as  discussed  previously,  it  is  not  possible  to  determine 
the  correct  positive  values  of  UT  approaching  separation.  It 
is  known  that  Ux  is  zero  at  xs.  It  is  the  purpose  of  the 
iteration  procedure  discussed  previously  to  find  the  correct 
location  of  xs. 

For  bodies  of  the  type  of  interest  herein,  the  imposition  of 
a conical  displacement  surface  corresponds  to  an  increasing  6* 
to  the  end  of  the  body  (figure  5c) . When  the  plume  or  the  sting 
surface  is  reached,  6*  begins  to  decrease.  Correspondingly,  UT , 
beginning  at  zero  at  xs  becomes  negative,  and  increases  in 
magnitude  until  the  end  of  the  body  is  reached.  For  a sting,  UT 
then  moves  back  toward  zero,  and  eventually  becomes  positive 
(reattachment).  At  some  point  after  reattachment,  the  calculation 
procedure  must  be  changed  back  to  the  ue  prescribed  mode  using 
equations  (47)  and  (48).  The  determination  of  the  point,  Xp  for 
this  switch  is  based  on  the  behavior  of  the  velocity,  ue,  produced 
by  the  boundary  layer  solution.  As  shown  in  figure  5,  ue  becomes 
nearly  constant  in  the  increasing  5*  region  between  xs  and 
the  end  of  the  body.  Then,  as  the  boundary  layer  proceeds  toward 
reattachment,  ue  decreases,  eventually  reaching  a minimum  a 
short  distance  downstream  of  reattachment.  If  <5*  is  prescribed 
beyond  the  point  where  the  minimum  velocity  is  reached,  a singu- 
larity develops  in  the  boundary  layer  solution.  To  avoid  this 
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singularity,  the  minimum  point  is  designated  as  Xp : and  the 
solution  is  switched  back  to  the  ue  prescribed  mode  (equations 
(47)  and  (48))  with  a ue  distribution  constructed  by  fairing 
the  viscous  velocity  from  the  minimum  value  into  the  inviscid 
distribution . 

The  procedure  just  described  was  found  to  produce  an  effec- 
tive body  shape  which  allowed  the  viscous  and  inviscid  velocity 
distributions  to  be  approximately  matched.  However,  to  achieve 
a matching  of  the  magnitude  of  the  two  velocities  and  locations 
of  the  minimum  velocity  on  the  sting,  it  was  found  to  be  necessary 
to  adjust  the  value  of  UT  at  Xp.  If  the  calculated  value 
of  Cf  at  that  point  is  less  than  0.001,  the  initial  value  of 
UT  for  the  calculation  downstream  of  Xp  is  adjusted  to  produce 
a value  of  Cf  of  0.001. 


6.0  EXHAUST  PLUME  ENTRAINMENT  MODEL 


The  flow  field  of  an  exhaust  jet  mixing  with  a!  subsonic 
external  flow  and  a separated  boattail  boundary  layer  is  shown 
in  figure  6.  Two  regions  can  be  distinguished  in  the  flow  field. 
In  region  1,  the  boundary  layer  merges  with  the  mixing  layer  at 
the  exhaust  plume  boundary  and  the  exhaust  plume  contains  a core 
of  inviscid  flow.  In  region  2,  the  turbulent  mixing  reaches  all 
the  way  to  the  axis  of  the  jet  so  there  is  no  inviscid  core  flow. 
In  this  report,  only  the  flow  in  region  1 is  considered  since 
for  the  kinds  of  flow  conditions  of  interest  the  entrainment  is 
expected  to  be  essentially  completed  in  the  length  of  the  in- 
viscid core. 

I 

6.1  ASSUMPTIONS 

The  analysis  is  based  on  the  following  assumptions: 

1.  The  flow  is  axisymmetric . 

2.  All  gases  obey  the  perfect  gas  law. 

3.  The  usual  boundary  layer  assumptions  are  applicable  to 
the  mixing  layer . 

4.  The  mixing  layer  is  fully  turbulent. 

5.  The  turbulent  Prandtl  and  Lewis  numbers  are  unity. 

6.  The  exhaust  jet  flow  is  supersonic  at  the  nozzle  exit. 

7.  Radial  pressure  gradients  are  negligible  in  the  mixing 
layer . 
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8.  Streamwise  turbulent  transport  mechanisms  are  negligible. 

9.  The  exhaust  jet  flow  is  isentropic. 

6.2  EQUATIONS  FOR  COMPRESSIBLE  TURBULENT  MIXING  LAYER 

The  mixing  layer  is  assumed  to  develop  about  a surface  de- 
fined by  the  radius  the  plume  would  have  if  it  was  completely 
inviscid  and  subject  to  the  pressure  distribution  calculated 
from  the  external  inviscid  flow.  The  basic  equations  are  assumed 
to  be  the  same  as  those  of  a turbulent  boundary  layer.  In  cylin- 
drical coordinates  for  axisymmetric  flow,  these  equations  have 
been  presented  previously  as  equations  (14)- (16).  For  the  en- 
trainment layer,  as  for  the  boundary  layer,  the  energy  equation 
(Eq.  16)  will  be  approximated  by  an  assumed  temperature-velocity 
relation, 

6.3  DEVELOPMENT  OF  INTEGRAL  METHOD  FOR  THE  MIXING  LAYER 


6.3.1  Integral  Equations 


For  the  entrainment  layer,  two  equations  are  derived  by 
eliminating  the  radial  velocity  terms  between  the  momentum  and 
continuity  equations  (14)  and  (15)  and  integrating  the  resulting 
equation  across  the  mixing  layer  in  two  strips.  No  compressibil- 
ity or  geometric  transformations  are  used.  Eliminating  the 
[v-ufr^y]  term  between  the  momentum  and  continuity  equations  and 
integrating  yields 


(pur)xdndr 


px  + rt$ur 


- P0r0Ivc-uc(rw)x](u-uc)  (67) 


where  it  has  been  assumed  that  the  shear  stress  term  vanishes  at 
rc , the  inner  boundary  of  the  mixing  layer.  Evaluating  equation 
(67)  at  r = rw  and  r = rP  yields  the  two  independent  integral 
equations 


f* 

| r [p  ^u-u^  ux  + 
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and 


f r[p(2u-ue)ux  + u<u-ue)px]dr  = 2peue(re"rc)  (u 


e>x 


- PGrc[vc"uc(rw)x] (ue_uc)  <69) 


where  the  pressure  has  been  eliminated  by  evaluating  equation  (15) 
at  r = re. 


The  external  boundary  of  the  mixing  region  is  defined  as 
follows : 


If 
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U1 

then 

re 

= 

rw 

(70) 


(71) 


where  <5U  is  the  thickness  corresponding  to  the  boundary  layer 
profile  at  the  end  of  the  boattail,  and  6C  is  the  thickness  of 
the  inner  part, of  the  mixing  layer. 

It  will  be  recalled  that  the  viscous-inviscid  interaction 
between  the  boundary  layer-mixing  layer  and  the  external  inviscid 
flow  requires  the  calculation  of  an  effective  boundary  shape  for 
the  inviscid  flow.  For  the  boundary  layer,  the  effective  shape 
is  defined  by  the  body  radius  plus  the  boundary  layer  displace- 
ment thickness,  equation  (55).  In  a similar  manner,  a displace- 
ment thickness  is  defined  for  the  exhaust  jet.  For  purely  in- 
viscid flow,  the  boundary  of  the  exhaust  jet  is  a streamline 
separating  the  nozzle  efflux  from  the  external  flow.  With  en- 
trainment and  turbulent  mixing,  this  line  is  no  longer  a stream- 
line. However,  in  order  to  produce  a continuous  and  smooth 
surface  for  the  external  inviscid  flow,  it  is  postulated  that  the 
effective  boundary  may  be  defined  by  a displacement  from  the 
inviscid  plume  boundary  in  the  same  manner  as  for  a solid  bound- 
ary. Thus,  an  additional  equation  is  derived  by  defining 

fIe 

5*  = I (l-pu/peue) (r/rw)dr  (72) 

r 

w 
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Differentiating  equation 
equation 


(r/rw)  <pu/peue)xdr  = - 


(72)  yields  a third  integral 

5x  • (rw)x(1'pwV'peV6Vrw1 


(73) 


6.3.2  Velocity  Profiles 

It  has  been  experimentally  observed  that  the  velocity  and 
temperature  profiles  in  jets  in  the  absence  of  entrainment  effects 
exhibit  self-similar  behaviors.  That  is,  appropriate  combinations 
of  the  radial  and  axial  coordinates  will  reduce  these  profiles 
to  functions  of  a single  variable  to  a high  degree  of  accuracy. 

In  the  present  case,  it  is  assumed  that  such  a condition  exists 
at  a large  distance  from  the  nozzle.  Near  the  nozzle,  the  flow 
profiles  are  assumed  to  be  composed  of  a linear  combination  of 
functions  representing  the  boundary  layer  profile  at  the  nozzle 
exit  and  the  equilibrium  jet  boundary  profile  far  downstream. 

Thus, 


u = (1-K)un  + Kug 


(74) 


where 


u^  is  a velocity  profile  proportional  to  the  boundary 

layer  profile  at  the  nozzle  exit  station  (both  internal 
and  external  boundary  layers) 

us  is  a velocity  profile  function  representing  the  down- 
stream profile 

and 

K is  a factor  to  be  determined. 


The  velocity  profiles  in  this  relation  are  defined  as  follows 


For  r < rw: 
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(75) 


where  y = r - r„ 
J w 
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For  r > r„ : 

w 


un  = lu'Vnue  (76> 

I 

where  (u/ue)  is  the  profile  of  the  external  boundary  layer  at 
the  end  of  tne  body,  A family  of  profiles  of  u for  several 
values  of  K and  ®c/rw  is  shown  for  illustration  in  figure  7. 
The  similarity  behavior  of  the  downstream  profile  is  assumed  to  ' 
be  approximated  by  a sine  function. 

Thus,  for  r - <5  < r < r.,  + 6 

w c w c 

us  = uc  + |<ue  - uc)[l  + sin(Tiy/26c)]  (77) 

and,  for  r > r,  + 

w c 


us  = % (78) 

With  these  definitions,  five  quantities  must  be  determined  to 
describe  the  flow.  Two  of  the  required  quantities,  rw  and  uc 
are  determined  by  the  approximate  model  of  the  flow  in  the  in- 
viscid  exhaust  jet.  The  remaining  three  quantities,  ue,  <$c  and 
K will  be  calculated  from  solutions  of  equations  (68)  , (69)  and 

(73)  combined  with  the  inviscid  external  flow  in  the  same  manner 
as  described  previously  for  boundary  layers. 

6.3.3  Temperature-Velocity  Relationship 

At  the  end  of  the  boattail,  the  temperature  parameter,  S is 
that  of  the  boundary  layer  as  described  previously.  For  the 
external  flow,  a function  Sn  is  thus  defined 

Sn  = [1  - (u/ue) ]SW  (79) 

Inside  the  nozzle,  a different  relation  is  needed  to  accommodate 
the  different  jet  total  temperature.  Assuming  the  linear  rela- 
tion between  the  total  temperature  and  the  velocity  also  exists 
inside  the  nozzle  boundary  layer,  the  following  relation  is  pro- 
posed for  r < rw; 


Sn  = [1  - {u/uc)]Sw  + Sc  (u/uc)  (80) 

For  fully  developed  plume  mixing  layers,  the  following 
temperature  profile  was  derived  by  Nielsen  et  al  (ref.  19) 
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(81) 


where 


c2  - 


m 


c 1 + m. 


(82) 


The  value  of  the  parameter  D has  the  following  effects; 

D = 0 gives  linear  T, -u  relation  and  quadratic  T-u 
relation. 

D = 1 gives  quadratic  Tfc-u  relation  and  linear  T-u 
relation. 


In  this  work,  D - 0 is  used  for  calculative  examples. 


In  order  to  express  all  terms  of  the  equations  in  terms  of 
a single  reference  quantity,  equation  (81)  is  rewritten  using 
the  external  flow  reference  conditions.  Thus,  the  similarity 
form  of  the  temperature  function  S becomes 


u -u  m 

S = S — — - + D •= — 

s c u -u_  1 + irv 


ru. 


u 


- 1 


ru 


<ue 


u. 


u. 


(83) 


For  the  region  immediately  downstream  of  the  nozzle  exit, 
a transitional  model  of  the  temperature  profiles  is  defined. 
Thus , 


S = (1  - K)Sn  + KSs  (84) 

The  density  profiles  are  related  to  the  velocity  profiles  through 
the  temperature  in  the  same  way  as  for  the  boundary  layer.  It 
is  assumed  that  the  mixing  in  the  entrainment  region  is  dominated 
by  the  external  flow,  since  for  separated  flow,  the  boundary 
layer  on  the  body  is  expected  to  be  much  thicker  than  the  boundary 
layer  inside  the  nozzle.  With  this  assumption,  the  ratio  of 
specific  heats  of  the  mixing  layer  is  taken  to  be  that  of  the 
external  flow.  The  density  profiles  are  then  obtained  from 

Pe/P  = (S+l)  (l+me)  - me(u/uj2  (85) 
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6.3.4  Turbulent  Shear  Stress 


Recalling  equations  (68)  and  (69) , it  is  noted  that  a re- 
lation is  required  for  the  shear  stress  along  the  inviscid 
reference  line,  r = rw.  In  equation  (68)  it  has  been  assumed 
that  the  shear  stress  is  related  to  the  radial  velocity  gradient 
through  an  eddy  viscosity  in  the  same  manner  as  for  the  boundary 
layer  on  the  body.  Also,  only  the  value  of  the  shear  at  the 
boundary,  r = rw  is  needed. 


As  for  the  velocity  and  temperature  profiles , the  eddy 
viscosity  for  a shear  layer  such  as  the  jet  boundary  far  down- 
stream of  the  nozzle  is  found  in  the  literature  to  have  a spe- 
cific form  different  from  that  of  the  boundary  layer.  In  the 
present  work,  the  eddy  viscosity  was  assumed  to  undergo  a trans- 
ition from  the  value  on  the  body  to  a value  determined  for  the 
downstream  flow  in  the  same  manner  as  for  the  velocity  and  temper- 
ature profiles,  namely. 


s = (l-K)v  + 2Ke 
o 


'f1+me] 

fU  +u  . % 

max  min 

l1+mJ 

^max  Umin-' 

0.9 


6 (u  -u  . ) 
c max  min 


(86) 


The  value  of  the  constant  e used  in  this  work  was 


e = 0.0009 


This  value  corresponds  to  the  correlation  of  Nielsen  et  al  (ref. 
19)  for  mixing  of  two  supersonic  streams. 

It  is  noted  that  the  shear  stress  term  in  equation  (68) 
contains  the  radial  velocity  gradient  as  well  as  the  eddy  vis- 
cosity. Differentiating  equation  (74)  , the  gradient  is 

Ur  = (1-K) (Un)r  + K(Us)r  (87) 


The  gradient  of  the  similarity  function  us  is  well  defined. 
However,  the  gradient  of  the  initial  boundary  layer  function  is 
not  well  defined.  The  boundary  layers  developing  on  the  boat- 
tail  surface  and  the  nozzle  interior  surface  cannot  be  expected 
to  have  a continuous  gradient  at  the  nozzle  exit.  However, 
immediately  downstream  of  the  nozzle,  a continuous  gradient  must 
exist  as  the  two  flows  begin  to  merge.  Since  the  jet  velocity 
is  expected  to  be  much  greater  in  magnitude  than  the  velocity  in 
the  recirculating  region  over  the  boattail , an  approximation  for 
the  velocity  gradient  at  the  inviscid  reference  line  was  derived 
by  assuming  that  the  velocity  near  the  boundary  r = rw  is 
describable  by  a sine  function  in  a similar  manner  to  the  down- 
stream similarity  function.  To  maintain  continuity  of  the 
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velocity  profile,  the  parameters  of  the  function  were  evaluated 
by  requiring  the  magnitude  and  gradient  of  the  velocity  profile 
to  be  continuous  at  rw-r/6c  = 0.1.  The  resulting  gradient  at 
r = rw  is 


(u  ) = -3 . 3931ttu  / 6 (B8) 

r rw  C c 

Calculations  have  indicated  that  this  relation  provides  a reason- 
able approximation.  The  solution  of  the  flow  equations  does  not 
appear  to  be  very  sensitive  to  errors  in  this  term. 


6.3.5  Equations  Solved 

The  integrands  of  equations  (68) , (69) , and  (73)  contain 

the  axial  derivatives  of  u and  p which  must  be  expressed  in 
terms  of  the  variables  ue,  6C,  and  K of  the  formulation.  The 
quantities  u and  p are  considered  to  be  functions  of  x and 
the  mixing  layer  variable  y = r - rw.  However,  the  axial  func- 
tional relationship  is  implicit  in  the  variation  of  the  other 
parameters  of  equations  (41) -(44)  so  that 


Then 


u (x,y)  = u(ue,6c,K,y) 
p(x,y)  = p{ue,6c,K,y) 


u„  = 


3u 

3u_ 


(ue'x 


3u 

35 


(5c>x 


3 u 
3K 


K. 


(89) 

(90) 


(91) 


3p 

3x 


(ue}x 


+ 


(6  ) + 
c x 


i£L  K 

3K  x 


(92) 


The  parameters  rw  and  uc  found  in  equations  (75) -(78)  are 
assumed  to  be  related  to  the  external  flow  velocity,  ue  and  do 
not  appear  explicitly  in  these  relations.  For  convenience,  the 
density  expression  is  rewritten  in  terms  of  the  density  ratio 
of  equation  (85).  Thus, 


Px  = <p/pe>  <Pe>x  - 


Pe (p/pe)  (pe/p)x 


(93) 


and 


(pe/c)x  = 


3u„ 


rf>. 


<ue>x 


35. 


&) 


<6c>x 


3K 
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where  it  can  easily  be  shown  that 


2p, 


«>e>x  - 


Y-l 


1+in 


1- 


1+m, 


iT  (1+me!<“e>x 
e 


(95) 


Due  to  the  iterative  nature  of  the  solution  procedure,  an 
exact  relation  cannot  be  derived  between  the  external  flow 
velocity,  ue  and  the  inviscid  plume  parameters,  rw  and  uc. 
However,  the  relationship  is  approximately  accounted  for  by  use 
of  the  assumption  that  the  solution  is  not  very  far  from  con- 
vergence so  that  both  the  external  and  the  plume  inviscid  flows 
are  subject  to  the  same  pressure  distribution,  pg.  The  external 
flow  is  isentropic  so  that 


1 


P u 
e e 


(96) 


Also,  for  an  inviscid  plume  which  is  nearly  uniform  and  axial, 
the  Bernoulli  equation  gives 


tc 


(97) 


so  that 


1 


pcuc 


(98) 


and 

3uc  _ 3uc/3Pe  = V^e 

3ue  " Sue/ape  pcuc 


(99) 


with  these  relations  and  the  definitions  of  equations  (75) -(78), 
the  required  derivatives  of  u can  be  obtained  as  follows 


3u 

= u -u 
3K  s n 


(100) 


7T-°  for  r > rw  + sc 
c 


when  < ^u 


(101) 


3u  3u 
= d-K)  ^ 

c c 


for  rw  - 3c  < r ‘ rw  + 6c 
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where 


3u  , . 

s = _ I (u-u) 

95c  4 


r-r. 


w 


cos 


1 jl) 

U aj 


and 


3u 

9? 


" 7 uc  (y/Sc)  (‘^c 


-6/7 


for  r„  - <5„<  r < r 
w c w 


= 0 for  r > r , and  6 <6 

w c u. 


= - y/6 


|u 

r 

u 

e 

iueJ 

for  r > r,,  and  6 >6 

w c u. 


J 


where  (u/ue)  1 denotes  differentiation  with  respect  to  y/<5( 
Continuing, 


3u 

3u_ 


3u 

(1-K)  (u/u6)  + ^ for  r > rw 


3u 

(1-K) (u/u  ) C 


3u 


o'  3u„  + K air  for  r < rw  j 


where 


3u 

■5—^  = 1 for  r > r + 6 
3u_  — w c 


= 3 [1  + sin(Try/26c)] 


3u  i 3u, 


1 - 


3u 


3u 


for  r -6  <r<r  +6 
w c — — w c 


Using  equation  (85)  in  equation  (94) 


(102) 


(103) 


(104) 


(105) 
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for  r > r, 
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and 


!i..s  u° : , P me  i 

3uc  C (U  -uj2  1+me  ue 
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(113) 


Also  , 
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as 

as 


as  au  as 

K sir  jr  + (1'K)  a<r 


(114) 


where 


as 

= 0 for  r > r and  6 < 5 

as  w c u. 
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,ZiL, 
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(115) 


Finally, 

51  = S - S 
3K  s n 

Equations  (68)  , (69)  , and  (73)  can  now  be  written,  after 

some  algebra,  in  the  form 


\ iA  ~r  n.  * j 

nl  x n2  c 


ing  definitions: 


w 


R _ o r_  _p_ 

3 p r p u 
e w *e  e 

_ 2 


P2  - ru  (/-)2(u  - uv) 


A - ( u ) = B 

n3  e x n 

(116) 

with  n = 1 , 2 or  3 

integrals  evaluated  from  the 

follow- 

> 

(117) 

,> 

(118) 

(119) 

V 

(120) 

uw) 

(121) 
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Thus , 


and 


where 


Also, 
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The  integrals  in  these  equations  are  evaluated  using  a 
21-point  trapezoidal  integration,  with  the  intervals  rc  to  rw 
and  rw  to  re  each  divided  into  10  equal  segments.  It  is  noted 
that  equations  (116)  are  similar  in  form  to  the  equations  solved 
for  the  boundary  layer  over  the  body.  Indeed,  a similar  approach 
is  taken  to  solving  the  entrainment  equations  as  was  used  for 
the  boundary  layer  equations.  At  the  end  of  the  body  (the 
beginning  of  the  entrainment  region) , the  value  of  6*  is  known, 
along  with  6U  and  ue.  A value  of  is  determined  from  a 

nozzle  boundary  layer  calculation.  The  value  of  K is  assumed  to 
be  zero.  The  variation  of  5*  is  specified  continuously  from  the 
-boundary  layer  into  the  entrainment  region,  while  ue,  <5C  and  K 
are  calculated  by  solving  equations  (116) . To  accomplish  this 
calculation,  the  first  two  of  equations  (116)  are  reduced  to  two 
equations  in  the  derivatives  of  <$cNand  ue  by  eliminating  the 
derivative  of  K with  the  third  equation.  The  value  of  K is 
determined  at  each  step  of  the  integration  of  the  solution  by 
solving  equation  (72)  directly.  When  the  point  of  minimum  ue 
is  reached  downstream  of  the  nozzle,  the  need  for  prescribing 
6*  is  eliminated.  From  that  point  to  the  end  of  the  calculation, 
ue  is  prescribed  as  for  the  boundary  layer  on  a sting  and  the 
values  of  K and  6C  are  calculated  from  the  first  two  of  equa- 
tions (116)  . 


7 . 0 RESULTS 

In  this  section,  the  theory  is  compared  with  experimental 
results  on  axisymmetric  boa ttail- sting  and  boattail-exhaust 
plume  configurations.  In  all  cases,  the  boundary- layer  was 
initialized  as  a laminar  boundary- layer  near  the  nose  of  the 
forebody  with  transition  specified  shortly  downstream  of  the 
initial  point.  Initial  boundary- layer  quantities  UT  and  <5^ 
were  obtained  from  a similarity  solution  transformed  to  conical 
coordinates  (see  section  10.4). 
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7.1  CIRCULAR  ARC  BOATTAILS  WITH  SOLID  PLUME  SIMULATORS 

Wind  tunnel  studies  of  the  configuration  shown  in  figure  8 
were  described  in  reference  20.  The  configuration  is  a cone- 
cylinder,  with  a circular  arc  boattail.  Two  of  the  boattails  are 
examined  here  to  illustrate  the  capabilities  of  the  calculative 
method.  The  first  example  is  for  a boattail  (configuration  2 in 
figure  8)  for  which  no  separation  is  evident  experimentally  at 
the  Mach  numbers  of  the  tests.  Surface  static  pressure  distribu- 
tions were  measured  along  the  boattail  and  along  the  sting  for 
free- stream  Mach  numbers  from  0.4  to  1.2.  Boundary- layer  transi- 
tion was  tripped  at  x/D  = .167.  Calculations  for  Mach  numbers  of 
0.6,  0.8,  and  0.9  were  initialized  as  described  previously.  On 
the  first  iteration  of  each  calculation  for  the  Mq  = 0.6  and  0.8 
cases,  a separation  point  was  indicated  on  the  boattail.  On 
subsequent  iterations,  the  separation  point  moved  aft  on  the 
boattail  until  it  eventually  moved  off  the  boattail  and  the  cal- 
culation subsequently  converged  as  an  attached  flow.  The  calcu- 
lation for  Mq  = 0.9  converged  without  separation.  Comparisons 
between  the  theory  and  measured  results  for  the  pressure  coeffi- 
cient distributions  at  Mach  numbers  of  0.6,  0.8,  and  0.9  are 
presented  in  figure  9.  The  theory  is  seen  to  be  in  good  agree- 
ment with  the  data  for  all  three  cases. 

It  will  be  recalled  from  section  5.6  that  the  initial  esti- 
mate of  the  separation  point  location  is  obtained  from  the  appli- 
cation of  the  plain- body  (the  body  with  no  boundary  layer) 
inviscid-flow  solution  to  the  boundary- layer  calculation,  with 
the  separation  point  taken  to  be  at  the  point  where  the  trans- 
formed shape  factor  reaches  a specified  value  in  the  adverse  ' 
pressure  gradient  on  the  boattail.  The  present  case  illustrates 
,that  the  calculative  method  also  has  the  capability 1 to  calculate 
attached  flows.  However,  for  such  flows  a separation  point  is 
often  found  on  the  first  iteration  because  of  the  steep  pressure 
gradients  associated  with  the  inviscid  flow  on  the  plain  body. 
Whether  the  subsequent  calculations  actually  proceed  to  move  the 
separation  point  off  the  boattail,  for  a case  in  which  that  is 
the  correct  result,  depends  somewhat  on  the  initial  value  of  xs, 
or  on  the  initial  value  of  H^.  If  ,the  initial  value  of  xs  is 
too  far  forward,  the  calculation  may  not  converge  at  all.  To 
move  the  initial  value  of  xs  back,  the  input  value  of  is 
increased.  For  the  cases  shown  in  figures  9a  and  9b,  a value  of 
H.^  of  1.3  was  appropriate,  while  for  the  case  shown  in  figure  9c, 
a value  of  1.5  was  required.  The  decision  as  to  whether  the 
solution  is  converged  is  based  primarily  on  the  magnitude  of  the 
rms  errors.  If  the  best  value  of  s produced  by  the  iteration  is 
significantly  greater  than  1.0  percent  (.01),  the  solution  can 
generally  be  improved  by  adjusting  H^ . 


40 


AEDC-TR-79-4 


The  second  example  is  for  a boattail  (configuration  1 in 
figure  8)  that  has  an  extensive  region  of  separated  flow  for  all 
Mach  numbers.  This  is  evidenced  by  the  pressure  distribution 
which  is  typical  of  separated  flows  on  such  configurations  and 
has  been  corroborated  by  oil  flow  studies  described  in  reference 
21.  Comparisons  between  the  theory  and  measured  results  for  the 
pressure  distributions  at  M0  = 0.6,  0.8,  and  0.9  are  presented 
in  figure  10.  For  the  two  lowest  Mach  numbers,  the  theory  is 
seen  to  be  in  excellent  agreement  with  the  data  on  the  boattail, 
with  a slight  discrepancy  on  the  solid  plume  simulator.  Also, 
the  calculated  location  of  the  separation  point  agrees  exactly 
with  the  oil  flow  determination  of  reference  21  for  the  Mq  =0.6 
case  of  figure  10a,  with  a small  error  noted  for  the  M0  = 0.8 
case  of  figure  10b.  For  the  M0  = 0.9  case  of  figure  10c,  the 
agreement  is  not  so  good.-  A fairly  large  discrepancy  is  noted 
in  the  pressure  coefficient  over  the  entire  boattail,  and  a 
larger  error  occurs  in  the  prediction  of  the  separation  location. 
For  that  case,  supersonic  flow  occurs  over  the  beginning  of  the 
boattail,  with  a shock  of  Mach  number  approximately  1.25  just 
ahead  of  the  separation  point.  Abeyounis  (ref.  21)  suggests  that 
for  this  configuration,  the  separation  for  free- stream  Mach  num- 
bers above  0.8  probably  occurs  at  the  shock.  Thus,  the  dis- 
crepancies noted  in  figure  10c  are  probably  due  to  the  fact  that 
the  present  theoretical  model  assumes  separation  occurs  behind 
the  shock.  Another  measure  of  the  accuracy  of  the  solution  is 
the  degree  of  convergence  of  the  calculated  viscous-inviscid 
interaction.  For  the  calculations  of  figures  10a  and  10b,  the 
iteration  converged  to  an  rms  error  of  1.0  percent,  which  has 
been  found  to  yield  generally  good  results,  while  results  shown 
in  figure  10c  had  an  rms  error  of  2,3  percent. 

As  for  the  unseparated  flows,  the  separated  flows  were  also 
found  to  be  somewhat  dependent  on  the  value  of  H-^  for  determining 
the  outcome  of  the  iterative  calculations.  In  the  two  subsonic 
cases  of  figures  10a  and  10b,  a value  of  H-^  of  1.15  was  used, 
while  for  the  transonic  case  of  figure  10c  a value  of  1-0  was 
used.  For  other  values  of  the  iteration  either  did  not  con- 

verge to  a sufficiently  small  rms  error,  or  developed  numerical 
difficulties  which  subsequently  stopped  the  calculations.  A 
value  of  Hi  = 1.3  is  usually  used  to  start  calculations  when 
nothing  is  known  of  the  solution.  Experience  indicates  that  if 
the  iteration  procedure  does  not  lead  to  a solution  with  an  rms 
error  of  1.0  percent  of  less,  it  may  be  possible  to  reduce  the 
error  by  adjusting  the  initial  value  of  Hi  and  starting  a new 
iteration  sequence. 
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Although  no  data  were  available  for  comparison,  other  bound- 
ary-layer quantities  of  interest  are  shown  in  figure  11.  The 
velocity  corresponding  to  the  pressure  of  figure  10b  is  shown  in 
figure  11a.  The  displacement  thickness,  also  shown  in  figure  11a, 
is  smooth  and  continuous,  with  an  abrupt  reversal  of  slope  at  the 
end  of  the  boattail,  consistent  with  the  prescribed  conical  sur- 
face of  the  effective  body  as  shown  in  figure  lib.  In  figure  lie,, 
the  skin-friction  coefficient  on  the  boattail  and  sting  is  shown. 
As  described  in  section  5.6,  the  skin  friction  is  negative  in  the 
region  between  xs  and  reattachment,  and  undergoes  an  adjustment 
at  Xp,  with  ue  and  6*  continuous. 

The  integrated  pressure  drag  on  the  boattails  is  shown  in 
figure  12.  The  calculated  values  at  Mq  = 0.6  and  0.8  are  in 
excellent  agreement  with  the  data  for  the  separated  flow  (figure 
12a) . For  higher  Mach  numbers,  the  agreement  is  not  so  good,  as 
would  be  expected. from  the  pressure  coefficient  comparisons  and 
the  previous  discussion.  Even  so,  the  calculations  for  Mq  = 0.9 
and  an  additional  calculation  at  Mq  = 0.96  indicate  the  method 
predicts  the  transonic  drag  rise  with  fair  accuracy. 

The  integrated  pressure  drag  on  the  unseparated  boattail  is 
shown  in  figure  12b.  The  error  in  the  calculated  drag  is  fairly 
large  for  the  Mo  = 0.6  case,  but  improves  considerably  for  the 
higher  Mach  numbers. 

7.2  CIRCULAR  ARC  BOATTAILS  WITH  COMPRESSED  AIR  PLUMES 

Some  of  the  wind  tunnel  tests  of  reference  20  were  conducted 
with  the  exhaust  plume  simulated  with  high  pressure  air.  Com- 
parison between  the  theory  and  data  for  the  pressure  coefficient 
distribution  on  the  boattail  of  configuration  1 (figure  8)  is 
presented  in  figure  13  for  exhaust  nozzle  total  pressure  ratios 
of  2.0  and  4.0  at  a free-stream  Mach  number  of  0.8  and  for  a 
nozzle  pressure  ratio  of  2.0  for  a Mach  number  of  0.6.  The 
agreement  between  the  theory  and  data  for  a pressure  ratio  of 
2.0  and  Mq  = 0.8  is  good  as  shown  in  figure  13a  even  though 
there  is  some  error  in  the  predicted  separation  point  location. 

The  rms  error  of  the  calculated  viscous  and  inviscous  solutions 
was  1.6  percent. 

The  comparison  for  a nozzle  pressure  ratio  of  4.0  and  Mq  = 

0.8  shown  in  figure  13b  indicates  good  agreement  for  the  pressure 
coefficient  distribution  and  somewhat  better  agreement  for  the 
separation  point  location  than  the  previous  case.  The  rms  error 
of  that  calculation  was  1.1  percent. 

The  comparison  for  a nozzle  pressure  ratio  of  2.0  and  Mq  = 

0.6  shown  in  figure  13c  indicates  fair  agreement  for  the  pressure 
coefficient  distribution  and  the  separation  point  location.  The 
rms  error  of  that  calculation  was  1.66  percent. 
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Plots  of  the  important  entrainment  region  parameters  for  the 
M0  = 0.8  case  are  presented  in  figure  14.  The  case  with  exhaust 
total  pressure  ratio  of  4.0  is  presented  in  figure  14a.  The 
plume  inviscid  reference  surface  represented  by  rw  is  seen  to 
expand  slightly  as  it  leaves  the  nozzle,  reaching  an  essentially 
cylindrical  shape  in  approximately  2 nozzle  diameters.  The  inner 
boundary  of  the  mixing  region,  represented  by  rc,  requires  more 
than  6 nozzle  diameters  to  reach  the  axis  of  the  plume.  The 
entrainment  parameter,  K increases  rapidly  at  first  and  then  more 
gradually  after  xp  is  reached  and  the  calculation  switches  from 
prescribed  6*  to  prescribed  ue.  Another  change  in  the  behavior 
of  K occurs  further  downstream,  associated  with  the  change  in  the 
definition  of  the  outer  edge  of  the  mixing  layer,  re,  as  defined 
by  equations  (70)  and  {71).  The  entrainment  reaches  completion 
(K  = 1.0)  asymptotically,  but  is  effectively  complete  in  less 
than  6 nozzle  diameters. 

The  entrainment  parameters  for  the  exhaust  pressure  ratio 
of  2.0  are  presented  in  figure  14b.  The  inviscid  plume  boundary, 
rw  is  cylindrical  for  the  entire  plume,  reflecting  the  fact  that 
the  plume  is  only  very  slightly  underexpanded.  Entrainment 
appears  to  proceed  at  a faster  rate  than  for  the  higher  pressure 
ratio  as  evidenced  by  the  fast  initial  rise  of  the  entrainment 
parameter,  K.  However,  the  entrainment  calculations  were  ter- 
minated after  approximately  2 nozzle  diameters  when  the  value  of 
K reached  a maximum  and  began  to  decrease.  The  reason  for  this 
behavior  is  not  fully  understood,  although  it  is  believed  to  be 
related  to  the  fact  that  the  jet  velocity  and  the  external  flow 
velocity  are  nearly  equal.  For  velocity  ratios,  uc/ue  close  to 
1.0,  it  is  expected  that  numerical  calculations  would  encounter 
difficulties  as  K approached  1.0.  For  K = 1.0  and  uc/ue  = 1.0, 
the  velocity  profiles  (Eqs.  (74)- (78))  become  constant  with 
radius,  and  the  equations  (68)  and  (69)  are  redundant.  The  ratio 
uc/ue  for  the  case  shown  in  figure  14b  is  approximately  1.2. 

Boattail  drag  coefficients  calculated  from  the  measured  and 
calculated  pressure  coefficient  distributions  of  figure  13  are 
shown  in  figure  15.  The  theoretical  drag  coefficients  for  the 
pressure  ratios  of  2,0  and  4.0  at  M0  = 0.8  are  in  excellent 
agreement  with  the  data  (figure  15a) . A third  case,  for  an 
exhaust  total  pressure  ratio  of  5.0  has  a calculated  drag  coeffi- 
cient slightly  higher  than  the  data.  The  rms  error  for  that 
case  was  3.2  percent.  The  theoretical  drag  coefficients  for  the 
Mq  = 0.6  case  shown  in  figure  15b  are  in  good  agreement  for  a 
pressure  ratio  of  4.0  but  the  agreement  for  lower  pressure  ratios 
is  poor.  The  rms  error  of  the  solutions  shown  in  figure  15b  was 
0.77  percent  for  Pt-/p0  = 4.0,  1.39  percent  for  pt_./p0  = 3-°  an<^ 

and  1.66  percent  for  pti/pQ  =2.0. 
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7.3  SUMMARY  OF  RESULTS 

The  calculative  method  has  been  found  to  be  applicable  to 
both  attached  and  separated  flows  on  boattails  with  solid  plume 
simulators  and  with  high  pressure  jets.  The  viscous-inviscid 
iteration  procedure  is  somewhat  sensitive  to  the  first  approxima- 
tion to  the  separation  point  location,  xs,  and  the  angle,  0S,  in 
that  the  iteration  will  not  converge  to  the  smallest  possible 
error  unless  the  initial  values  are  close  to  the  correct  values. 
This  does  not  appear  to  be  a major  limitation,  since  adequate 
starting  conditions  can  usually  be  determined  by  trial  and  error 
in  a few  trials.  Comparisons  with  experimental  data  indicate 
that  the  entrainment  model  and  the  viscous-inviscid  iteration 
procedure  provide  an  accurate  engineering  method  for  predicting 
boattail  flow  fields  for  moderately  underexpanded  exhaust  flows, 
and  for  boattails  with  solid  exhaust  plume  simulators. 


8.0  COMPUTER  PROGRAM  ORGANIZATION 


In  this  section,  the  general  organization  of  the  programs 
will  be  described.  Specific  information  on  data  required  for 
input  and  data  developed  for  output  will  be  described  in  sections 
9,  10,  and  11.  Sample  Job  Control  Card  decks  for  a typical  IBM 
370  installation  are  presented  in  section  12. 

The  overall  program  consists  of  a mainline  program  {program 
1)  and  three  main  subprograms  each  consisting  of  several  sub- 
routines. The  first  main  subprogram  is  the  inviscid-f low  pro- 
gram (program  2) . It  is  a modified  version  of  the  program 
described  in  reference  3.  The  second  main  subprogram  is  the 
boundary- layer- entrainment-layer  program  (program  3} . It  is 
based  on  the  integral  theory  described  previously.  The  third 
subprogram  is  the  inviscid  exhaust  plume  program  (program  4) . 

The  mainline  program  controls  the  iteration  between  the  other 
three  programs.  Either  of  the  first  two  subprograms  may  be  used 
separately,  without  iterating  by  appropriate  choice  of  the  input 
parameters. 

Both  the  viscous-flow  program  and  the  inviscid-f low  programs 
require  some  punched  card  input  and  may  require  some  input  data 
from  disc  or  tape  data  files.  The  data  files  must  be  identified 
by  specific  Logical  Unit  numbers.  Each  program  in  turn  produces 
new  data  files  and  printed  output.  The  general  relationship  of 
the  programs  and  the  various  data  files  are  shown  in  figure  8. 

The  specific  Logical  Unit  numbers  required  for  input  and  output 
are  listed  in  Table  I.  Two  Logical  Unit  numbers  are  associated 
with  each  data  file  shown  in  figure  16.  One  unit  is  used  for 
input,  the  other  for  output. 
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Program  2 requires  initially  data  from  cards  describing  the 
free-stream  conditions,  the  computational  mesh  and  the  body  shape. 
Alternately,  the  program  can  accept  the  input  body  shape  from 
data  file  1.  It  can  also  accept  an  initial  solution  for  the 
perturbation  velocity  potential  from  another  data  file  (data 
file  2) . Program  2 produces  printed  output  lists  of  the  appro- 
priate flow  field  quantities,  quantities  describing  the  config- 
urations and  the  computational  mesh  and  several  data  files.  Data 
file  2 is  rewritten  using  the  new  solution  for  the  potential.  A 
third  file  (data  file  3)  is  written  containing  the  distribution 
of  the  axial  velocity  component  at  the  inviscid  boundary  for  use 
by  the  viscous  flow  program. 

Program  3 requires  initially  the  free-stream  conditions  and 
gas  constants  as  well  as  parameters  describing  the  shape  of  the 
surface  over  which  the  boundary  layer  is  flowing.  On  the  first 
iteration  of  a viscous-inviscid  interaction,  the  surface  shape 
is  the  same  as  that  for  program  2.  On  subsequent  iterations, 
the  body  shape  for  the  boundary-layer  program  remains  the  same, 
except  when  it  is  modified  by  the  calculation  of  a new  plume 
shape,  while  that  for  the  inviscid  flow  is  modified  by  the 
addition  of  the  boundary  layer.  Program  3 can  also  accept  data 
from  data  files  as  optional  input.  The  boundary- layer- edge 
velocity  distribution,  ue , can  be  input  from  data  file  3,  pro- 
duced by  program  2.  The  distribution  of  the  body  shape  augmented 
by  the  displacement  thickness  can  be  input  from  data  file  4 . 

That  file  differs  from  data  file  1 because  it  contains  the  raw 
data  for  &*  + rw  versus  x as  calculated  by  the  boundary-layer 
program  while  data  file  1 contains  the  shape  adjusted  according 
to  equation  (4)  and  interpolated  to  the  x stations  of  the  original 
input  shape.  If  program  3 is  being  used  separately  from  program 
2 (i.e.,  without  iterating),  two  additional  options  are  available 
for  card  input.  Either  the  boundary-layer-edge  velocity,  ue,  can 
be  input  as  mentioned  previously,  or  the  displacement  thickness, 
6*,  may  be  input.  These  options  are  described  more  fully  in 
sections  10.1  and  10.2. 

Program  3 produces  as  output  lists  of  the  boundary  layer  and 
flow  quantities  as  they  are  calculated  along  the  body.  In 
addition,  program  3 produces  an  updated  version  of  data  file  4 
and  the  augmented  body  shape  (data  file  1)  required  by  program  2. 
Data  file  4 also  contains  a list  of  the  velocity  ratio  ue/ue 
corresponding  to  the  boundary  layer.  ° 

Program  3 can  be  used  in  a two-dimensional  boundary  layer 
mode  if  desired.  However,  calculation  of  a viscous-inviscid 
interaction  or  exhaust  plume  entrainment  can  only  be  done  for 
axisymmetric  cases  with  the  present  inviscid  program. 
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All  card  input  pertaining  to  programs  2 and  3 is  input 
through  an  input  subroutine  called  by  program  1.  Program  1 also 
produces  a file  of  the  quantities  needed  to  restart  and  continue 
the  calculation  if  the  calculation  should  terminate,  before  all 
iterations  are  completed.  These  data  are  stored  on  data  file  5. 
Detailed  instructions  regarding  restarting  are  presented  in 
sections  12.2.3  and  13.2. 


9.0  INPUT  TO  THE  PROGRAMS 


The  data  required  by  the  programs  generally  fall  into  three 
categories:  (1)  geometrical  data;  (2)  flow  field  data;  and  (3) 

control  parameters.  The  control  parameters  are  indices  for 
specifying  options  and  iteration  counters.  It  will  be  noted  by 
comparison  with  reference  3 that  a number  of  input  quantities 
required  for  the  inviscid-f low  program  have  been  eliminated  in 
the  present  version.  This  has  been  done  by  incorporating  the 
calculations  required  to  obtain  some  of  the  quantities  into  the 
present  code  or  by  simply  defining  fixed  values  which  have  been 
found  to  be  successful.  Specifically,  a value  of  1.4  is  used 
for  the  initial  value  of  the  subsonic  relaxation  factor,  a value 
of  0.1  is  used  for  the  initial  value  of  the  supersonic  relaxation 
factor  and  a value  of  1.3  is  used  for  the  exponent  in  the  normal 
coordinate  stretching  function.  Also,  it  is  assumed  that  the 
computational  grid  has  equal  step  sizes  in  both  coordinate  direc- 
tions at  the  nose  of  the  body.  Finally,  the  number  of  relaxation 
steps  allowed  in  the  inviscid  program  is  fixed  at  20  for  the 
first  four  interaction  iterations  and  is  changed  to  40  and  80 
after  four  and  eight  iterations,  respectively. 

The  general  requirement  of  the  input  data  is  that  the  tabu- 
lar lists  of  the  various  distributions  required  represent  smooth 
curves.  This  is  especially  true  of  the  list  of  body  shape  coor- 
dinates. The  inviscid  program  uses  cubic  splines  to  fit  the 
input  coordinates,  so  those  coordinates  must  accurately  represent 
a smooth  curve  with  continuous  second  derivatives. 

9 . 1 TABULAR  FORM 

The  input  data  required  for  calculating  transonic  viscous, 
inviscid-f low  interactions  consist  of  several  punched  cards  con- 
taining parameters  describing  the  free-stream  flow  conditions, 
the  computational  mesh  for  the  inviscid  calculation,  initial 
values  for  the  viscous  flow  calculation,  and  certain  options 
that  are  available  in  the  program.  A dictionary  of  the  input 
data  is  presented  in  the  next  section.  Table  II  shows  the  input 
variables  as  they  are  to  be  punched  on  the  data  cards.  More 
detailed  explanation  of  the  requirements  for  the  inviscid-f low 
program  (program  2)  are  presented  in  ref.  3 and  are 'not  repeated 
herein. 
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9.2  DICTIONARY  OF  INPUT  VARIABLES 

The  variables  required  for  input  on  punched  cards  are 
defined  in  this  section  in  the  order  in  which  they  are  required. 
Additional  details  on  the  format  of  the  punched  data  are  given 
in  Table  II.  The  first  three  cards  of  any  input  data  deck  con- 
tain a description  of  the  case  being  calculated.  Any  or  all  of 
these  three  cards  may  be  blank,  but  all  three  are  required.  The 
remaining  variables  in  Table  II  are  as  follows: 

NRSTRT  Integer  indicating  whether  calculation  is  being  restarted 
to  continue  a previous  calculation.  Only  valid  for 
interaction  calculations  (LPROG  = 0.  See  below) 

= 0 Start  from  zero.  Input  all  quantities  on  cards 
or  data  files  as  required. 

> 0 Restart.  Input  data  file  5 (Logical  Unit  11) 

containing  data  from  previous  iteration  plus  all 
other  input  data  files  plus  other  card  data  as 
indicated  in  Table  II. 

NPRINT  Integer  indicating  quantity  of  output  to  be  printed 
(see  section  11.1). 

= 0 Minimum  output. 

= 1 All  output. 

N3  Integer  indicating  whether  restart  is  to  begin  with  new 

values  of  xs  and  9S  (see  section  13.2). 

ILIM  Integer  number  of  seconds  corresponding  to  estimated 

length  of  calculation.  When  the  time  from  the  beginning 
of  the  calculation  is  within  30  seconds  of  this  limit, 
the  final  restart  files  will  be  written.  This  assures 
that  the  job  will  not  terminate  while  writing  such  files. 

I I TER  Integer  determining  size  of  perturbation  of  xg  and  0S. 

The  initial  step  sizes  used  will  be  those  given  by 
equations  (8)  and  (9)  divided  by  IITER. 

LPROG  Integer  indicating  level  of  calculation. 

= -1  Inviscid  flow  only. 

= 0 Viscous-Inviscid  interaction. 

= 1 Boundary  layer  only. 


47 


A E DC-T  R-79-4 


N1 

N2 

IBL 

IUNIT 


MIT 

GAM 

AMINF 

IXY 

XO,YO 


Integer  iteration  counter  for  inner  viscous- inviscid 
iteration. 

Integer  iteration  counter  for  xg-0s  iteration  presently 
limited  to  a maximum  of  20. 

Integer  indicating  how  interaction  calculations  are  to 
begin. 

= 3 Start  with  inviscid  flow. 

= 0 Start  with  boundary  layer. 

Integer  indicating  which  value  of  the  gas  constant,  RGAS , 
and  the  constants  in  Sutherlands  temperature-viscosity 
relation  are  to  be  used.  The  choice  depends  on  whether 
air  is  the  gas  being  calculated  and  the  units  of  the 
input  quantities. 

= 1 Input  units  must  be  pounds,  feet,  seconds,  and 
°R. 


= 2 Input  units  must  be  pounds,  inches,  seconds,  and 
°R. 

= 3 Input  units  must  be  newtons,  meters,  seconds,  and 

°K. 


For  another  gas,  or  other  units,  put  in  anything  for 
IUNIT  and  put  in  nonzero  values  of  VISC,  RGAS,  and  SCON. 
Any  units  are  allowed.  The  basic  rule  is  that  all  input 
quantities  must  be  consistent  with  regard  to  units. 

Number  of  relaxation  cycles  allowed  for  the  inviscid-f low 
program  {Program  2) . If  a value  of  zero  is  input,  a 
value  of  20  is  used.  Otherwise,  the  input  value  is  used. 

Ratio  of  specific  heats  of  the  external  inviscid  flow. 

Free- stream  Mach  number. 

Integer  number  of  values  of  coordinate  pairs,  XO,YO,  to 
be  input  for  inviscid  body  shape.  If  IXY  = 0,  the 
required  shape  must  be  input  from  data  file  1 (Logical 
Unit  14) . Maximum  value  is  200. 

Axial  and  radial  coordinates  of  body  shape  for  inviscid- 
flow  calculation,  2 per  card. 
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The  next  series  of  variables,  items  7 and  8 in  Table  II,  are  for 

the  inviscid-f low  program  (Program  2) . Detailed  information  on 

how  these  data  are  to  be  obtained  is  contained  in  reference  1. 

IMAX  Number  of  grid  lines  in  the  tangential  direction;  1=1 
is  the  forward  stagnation  line,  I = IMAX  is  the  rear 
stagnation  line  for  closed  bodies  and  downstream  infinity 
for  open  bodies.  For  each  grid  refinement  IMAX  is 
increased  such  that  IMAXjjg^  = 2(IMAXqld)  “ 1*  The  pres- 
ent limit  on  IMAX  is  81.  Instructions  for  changing  this 
limit  appear  as  comments  in  the  program  listing  (sub- 
routine ONEO) . The  grid  refinement  option  is  generally 
not  used  for  viscous-inviscid  interaction  calculations. 

JMAX  Number  of  grid  lines  in  the  normal  direction;  J = 1 

corresponds  to  an  infinite  distance  from  the  body  and 
J = JMAX  is  on  the  body.  The  same  formula  and  limit  that 
apply  to  IMAX  also  apply  to  JMAX. 

MHALF  Number  of  grid  refinements  to  be  done.  For  interaction 
calculations  a value  of  zero  should  be  used  with  a grid 
fine  enough  for  adequate  resolution. 

KLOSE  Body  type. 

= 0 Open  body  ’(i.e.,  one  with  a sting  or  wake). 

= 1 Closed  body. 

LREADP  Integer  indicating  whether  initial  estimate  of  potential 
distribution  is  to  be  input  from  data  file  2 (Logical 
Unit  13) . 

= 0 No. 

= 1 Yes. 

If  any  one,  or  all  of  the  next  four  input  quantities  are  input  • 

as  zero,  the  pre-programmed  values  are  used. 

DNDZO  Step  size  of  the  normal  coordinate  at  the  body.  The  pre- 
programmed value  is  7 percent  of  the  maximum  body 
diameter . 

XIXM  Value  of  the  computational  coordinate,  X,  at  the  match- 
ing point  of  the  two  stretching  functions  used  in  the 
finite-difference  scheme  (see  ref.  3) , for  open  bodies 
only.  Since  X varies  from  zero  to  one,  XIXM  is  the 
fraction  of  the  total  number  of  grid  points  which  will 
be  in  the  first  stretching  region  (ahead  of  x^) . The 
pre-programmed  value  is  0.75. 
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XM  Axial  location,  xm,  (in  physical  coordinates)  of  the 

junction  (or  matching  point)  between  the  two  tangential 
stretching  functions,  for  open  bodies  only,  see  refer- 
ence 3.  Must  be  less  than  XO(IXY) . This  parameter  is 
used  to  concentrate  computational  mesh  points  in  a cer- 
tain region.  The  usual  approach  for  interaction  calcu- 
lations (the  pre-programmed  value)  is  to  let  xm  be  equal 
to  the  length  of  the  body  to  the  beginning  of  the  sting, 
or  plume , XBT . > 

DSDXIM  Step  size  of  the  tangential  coordinate  at  the  junction 
between  the  two  tangential  stretching  functions.  The 
pre-programmed  value  is  8 percent  of  the  length  of  the 
afterbody  (XBT-XZNEW) . 

XBT  Length  of  the  body. 

DMAX  Maximum  body  diameter. 

XZNEW  Length  of  forebody.  Also,  this  is  the  axial  location  at 
which  boundary- layer  calculations  will  begin  after  four 
iterations  have  been  calculated.  The  usual  procedure  is 
to  start  the  boundary- layer  calculations  close  to  the 
nose  of  a long  body  at  XZ  (see  item  17  in  Table  II)  and 
then  after  four  iterations  move  the  starting  point  to 
XZNEW  (XZNEW  > XZ) . In  subsequent  iterations,  the  bound- 
ary layer  does  not  change  for  X < XZNEW.  For  long 
slender  bodies  with  boattails,  XZNEW  should  be  the 
beginning  of  the  boattail. 

PLUMIN  Logical  variable  indicating  whether  an  exhaust  plume  is 
being  calculated. 

GAMAP  Ratio  of  specific  heats  of  exhaust  gas. 

PTPPFS  Ratio  of  nozzle  total  pressure  to  free  stream  static 
pressure. 

AMP  Nozzle  exit  Mach  number. 

THETAP  Nozzle  exit  divergence  angle. 

TTP  Nozzle  total  temperature. 

GMP  Ratio  of  gas  constant  of  air  to  that  of  the  exhaust  gas. 

DELCZ  Thickness  of  nozzle  internal  boundary  layer  at  the 

nozzle  exit. 
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The  remaining  input  variables  are  related  to  the  boundary- layer 
program  alone  or  to  the  viscous-inviscid  interaction  method. 

LSEP  Logical  variable  indicating  whether  the  location  of 
separation  (XSEP)  is  known  (TRUE)  or  not  (FALSE) . 

XSEP  Axial  location  of  separation.  Put  in  only  if  LSEP  = TRUE. 

THETS  Angle  of  £*  surface  with  the  boattail  tangent  at  xs. 

Put  in  only  if  LSEP  = TRUE.  A value  of  0.0  will  cause 
the  program  to  calculate  a value  using  equation  (3) . 

IOPT  Integer  indicating  the  mode  of  the  calculation. 

=1  ug  is  to  be  input. 

=2  6*  is  to  be  input. 

Put  in  a value  of  1 for  starting  an  interaction  calcu- 
lation. 

K Integer  indicating  whether  flow  is  axisymmetric . 

= 0 Two  dimensional . 

= 1 Axisymmetric. 

LVARl  Integer  indicator  for  method  of  input  of  u=  when 
IOPT  *»  1. 

= 0 Input  u (dimensional)  on  cards. 

= 1 Input  ue/u@o  on  cards . 

= 2 Input  u„/u„  from  data  file  3 on  Logical  Unit  12. 

eo 

LSHAPE  Integer  indicating  option  for  calculating  all  initial 
conditions  (I.C.)  except  ug  (see  section  10.4). 

= 0 Input  initial  values  per  LIC. 

= 1 Calculate  I.C.  for  flat  plate. 

= 2 Calculate  I.C.  for  cylinder. 

= 3 Calculate  I.C.  for  cone. 

LIC  Integer  indicating  initial  condition  options  for  IOPT  = 1 

and  LSHAPE  =0. 
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= 1 Put  in  CFCl  and  DELTA 1 . 

= 2 Put  in  CFCl  and  DELSTl . 

LDSTAR  Integer  indicating  whether  a file  of  6*  + r is  to  be 
input. 

=0  No  input. 

= 1 File  of  6*  + rw  versus  x is  required  on  Logical 
Unit  15  (data  file  4). 

LSHPBL  Integer  indicating  whether  body  shape  is  to  be  input  for 
boundary- layer  calculations. 

= 0 XRL  and  RL  are  assumed  to  be  the  same  as  XO  and 
YO.  This  is  usually  the  case  when  starting  an 
interaction  calculation. 

= 1 XRL  and  RL  will  be  required. 


NOTE:  The  next  three  variables,  NVAR,  XVAR,  VAR,  are  only 

required  on  cards  if  LPROG  = 1 and  LVARl  = 0 or  1. 


NVAR  Integer  indicating  the  number  of  values  to  be  input  for 
the  prescribed  variable  (u  or  5*) . Maximum  value  is 
100. 

XVAR,  Axial  location  and  value  of  prescribed  variable  as 
VAR  follows: 


IOPT  = 1 and  LVARl  = 0,  VAR  = u 


EL 

PT 

TT 

TWONTT 

vise 


IOPT  = 1 and  LVARl  = 1,  VAR  = u /u 
IOPT  = 2,  VAR  =6* 


'o 


Reference  length.  Needed  if  input  data  lengths  are  non- 
dimens ionali  zed.  If  lengths  are  dimensional,  put  in 
EL  = 1.0. 

2 2 2 

Total  pressure,  pt  (lb/ft  ),  (lb/in  ),  or  (newton/m  ). 
Total  temperature,  T^  (°R)  or  (°K) . 


Ratio  of  body  surface  temperature  to  total  temperature, 

VTf 

Constant  \ in  Sutherlands  formula  for  viscosity. 
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RGAS 


SCON 


DFACT 

XZ 

RLEN 


T3/2 

U ~ A T + T 

S 

If  one  of  the  programmed  values  is  acceptable/  put  in  a 
value  of  0.0.  The  value  used  will  then  be  determined  by 
the  value  of  IUNIT  on  item  number  3 as  follows: 

I UNIT  VISC 

1 2.27(10~8)  lb  sec/ft2  {°R)1/2 

2 1.5764  (IQ-10)  lb  sec/in2  (°R)1/2 

3 1.4582  (10-6)  Newton  sec/m2  (“K)^2 

Gas  constant.  If  one  of  the  programmed  values  is 
acceptable,  put  in  a value  of  0.0.  The  value  used  will 
be  determined  by  the  value  of  IUNIT  as  follows: 

IUNIT  RGAS 

1 1716.0  ft2/sec2  °R 

2 247104.0  in2/sec2  °R 

3 286.96  m2/sec2  °K 

Constant  Ts,  in  Sutherlands  viscosity  law  (see  definition 
for  VISC  and  equation  31) . If  one  of  the  programmed 
values  is  acceptable,  put  in  a value  of  0.0.  The  value 
used  will  be  determined  by  the  value  of  IUNIT  as  follows. 

IUNIT 

1 198.6  °R 

2 198.6  °R 

3 110.333  °K 

Relaxation  factor  for  adding  6*  to  body,  the  factor  a 
in  equation  (4).  Usual  value  is  0.5. . If  a value  of 
0.0  is  put  in,  a value  of  0.5  will  be  used. 

Axial  location  of  beginning(  of  boundary- layer  calculation. 

Axial  location  of  end  of  boundary- layer  calculation. 
Usually  at  least  one  maximum  diameter  larger  than  XBT. 

For  a plume  RLEN  is  automatically  extended  to  match  the 
inviscid  grid. 
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XT 

DXP 

HLIM 

CFC1 

t 

DELTA1 

DELST1 

UE1 

DUEDX 

NR 

XRP , RL 


Axial  location  of  transition  from  laminar  to  turbulent 
boundary- layer  flow. 

Axial  interval  at  which  velocity  profiles  are  to  be 
printed.  If  a value  of  0.0  is  put  in,  no  profiles  are 
printed.  In  any  case,  profiles  will  not  be  printed 
more  often  than  the  output  step  size  for  individual 
output  quantities  (see  section  11.1). 

Limit  value  of  to  indicate  separation  as  discussed 
in  section  5.6.  Input  of  a blank  or  a value  of  0.0  will 
cause  a value  of  1.3  to  be  used.  To  move  the  initial 
value  of  xs  forward,  decrease  HLIM.  To  move  xs  down- 
stream, increase  HLIM.  A value  greater  than  4.0  will 
allow  the  calculation  to  proceed  to  a separation  singu- 
larity if  one  should  occur.  The  calculation  will  ter- 
minate at  that  point.  See  section  7 for  a discussion 
of  the  use  of  HLIM  in  sample  calculations. 

Value  of  skin-friction  coefficient  at  initial  boundary- 
layer  station  (see  section  10.4). 

Value  of  boundary- layer  thickness,  6 (compressible) , at 
initial  boundary- layer  station  (see  section  10.4). 

Value  of  boundary- layer  displacement  thickness,  6*,  at 
initial  boundary- layer  station  (see  section  10.4). 

Value  of  boundary-layer-edge  velocity  ue  at  initial 
boundary-layer  station. 

Value  of  boundary- layer-edge  velocity  gradient  at  initial 
boundary- layer  station. 

Integer  number  of  values  of  XRP  and  RL  to  be  input  for 
body  shape.  If  NSHPBL  = 0,  this  is  assumed  to  be  the 
same  as  IXY.  Maximum  value  is  200. 

Axial  and  radial  coordinates  of  body  shape  for  boundary- 
layer  calculation.  If  LSHPBL  = 0,  these  are  assumed 
to  be  the  same  as  XO  and  YO,  respectively.  For  two- 
dimensional  configurations,  these  represent  the  x and  y 
coordinates  of  a surface  measured  from  a reference  plane. 
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10.0  PROGRAM  OPTIONS 


Several  optional  inodes  of  calculation  are  available  through 
the  input  parameters.  A description  of  the  options  and  the 
corresponding  values  of  the  pertinent  parameters  follows. 

10.1  BOUNDARY- LAYER  OPTION 

To  use  only  the  boundary- layer  program,  put  in  the  three 
card  description,  then  all  values  on  the  fourth  input  data  card. 
Item  2 in  Table  II,  should  be  zero.  Then  the  first  value  on  the 
fifth  input  data  card  should  be: 

LPROG=l  (see  Table  II,  Item  3) 

Of  the  remaining  variables  shown  on  item  number  3 in  Table 
II,  only  I UNIT  is  required.  The  remaining  cards  would  be  those 
corresponding  to  item  4 and  items  13  to  21  as  described  in  Table 
II.  With  this  option,  the  user  has  the  choice  of  specifying 
either  the  free-stream  velocity,  ue,  or  the  displacement  thick- 
ness, S*,  through  the  variables  NVAR,  XVAR,  and  VAR  on  items  14 
and  15.  The  boundary-layer  calculation  can  be  restarted  at  any 
station  by  inputting  the  values  of  the  variables  at  that  station 
as  listed  in  the  output.  The  calculated  list  of  6*  + rw  and 
ue/ue  will  be  written  on  Logical  Unit  10  (data  file  4)  when  the 

calculation  terminates. 

Another  method  is  also  available  for  calculating  the  bound- 
ary layer  alone.  The  boundary- layer  step  of  a viscid- inviscid 
iteration  can  be  executed  separately.  The  appropriate  values  on 
the  fifth  card  would  be: 

LPROG  = 0 
Nl  = 1 
N2  = 21 
IBL  = 0 

IUNIT  = 1,  2,  or  3 

This  option  requires  input. of  all  quantities  as  though  the 
iterative  sequence  were  to  be  completed.  With  the  values  just 
described,  only  the  boundary  layer  will  be  calculated  and  then 
the  run  will  be  terminated.  If  it  is  desired  to  continue  the 
iteration,  simply  make  N2  less  than  21.  The  free-stream  velocity 
distribution  must  be  provided  on  Logical  Unit  12  for  this  case. 
All  other  optional  inputs  are  the  user's  choice. 
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10.2  INVISCID-FLOW  OPTION 

To  use  only  the  inviscid-f low  program,  put  in  the  three  card 
description,  then  all  values  on  the  fourth  card  should  be  zero. 
Then  put  in  LPROG  = -1  on  the  fifth  card  (item  number  3 in  Table 
II) . Of  the  remaining  values  on  that  card,  only  IUNIT  is 
required  and  that  quantity  is  required  only  if  PLUMIN  is  TRUE. 
After  the  first  five  cards  only  the  data  for  items  4 to  9 or  10 
as  described  in  Table  II  and  section  9.2  are  required  for  this 
option. 

10.3  VISCID-INVISCID  ITERATION  OPTION 

To  use  both  the  boundary- layer  and  the  inviscid-f low  programs 
iteratively,  put  in  LPROG  = 0 and  all  other  quantities  as  appro- 
priate. Such  iterations  can  be  started  with  only  the  body  shape 
and  free- stream  flow  quantities  known  and  may  be  restarted  to 
continue  a prematurely  terminated  iteration.  Several  options  are 
available  to  the  user  for  restarting  an  unfinished  iteration. 

See  section  13.2  for  an  example  of  restarting.  The  simplest 
option  is  to  put  in  NRSTRT  = 1 as  the  first  value  on  the  fourth 
input  data  card  and  to  provide  the  required  input  data  files  on 
Logical  Units  11,  12,  13,  14,  and  15  (see  fig.  16  and  Table  I). 

The  only  other  data  required  for  restarting  are  the  three-card 
description  of  the  case.  The  calculation  then  proceeds  from 
where  the  previous  iteration  stopped.  Another  method  of  restart- 
ing would  be  to  omit  the  restart  file  and  put  in  NRSTRT  * 0.  The 
user  can  then  vary  any  of  the  other  input  quantities,  using  the 
data  files  or  punched  cards  as  desired.  Note  that  the  calcula- 
tion terminates  when  N2  reaches  a value  of  21.  The  value  of  Nl 
increases  continuously  throughout  the  calculation  while  the  value 
of  N2  is  reset  to  1 each  time  the  xs-0s  cycle  finds  a minimum 
error  without  satisfying  the  convergence  criterion. 

10.4  BOUNDARY- LAYER  INITIAL  CONDITIONS 

Initial  values  of  boundary- layer  quantities  can  be  obtained 
in  several  ways.  The  user  can  obtain  values  of  the  skin-friction 
coefficient,  Cf,  and  either  the  boundary- layer  thickness,  6,  or 
the  displacement  thickness,  6*.  These  are  shown  on  item  19  in 
Table  II.  The  appropriate  values  LSHAPE  = 0 and  LIC  = 1 or  2 
are  punched  on  the  card  corresponding  to  item  13.  For  the  case 
when  no  other  source  of  this  information  is  available,  formulas 
have  been  included  in  the  program  based  on  the  Blasius  solution 
for  laminar  boundary  layers  and  based  on  one- seventh  power  law 
velocity  profiles  for  turbulent  flows.  These  formulas  are  only 
available  if  ue  is  being  specified  {IOPT  = 1) . The  basic 
formulas  calculate  Cf  and  6 in  the  transform  plane  (incompress- 
ible, two-dimensional) . The  formulas  are  as  follows: 
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Laminar  Flow 


6 . 
i 


0.664 


Turbulent  Flow 


C 


f . 
i 


0.0592 


0.37  x 


'1-0.2 


x 


(127) 


(128) 


(129) 


(130) 


These  formulas  provide  initial  values  for  boundary  layers  on  flat 
plates.  They  are  chosen  by  inputting  LSHAPE  = 1:  For  other 
geometries,  the  value  of  the  x coordinate  is  transformed.  Thus, 
LSHAPE  = 2 chooses  the  values  for  a circular  cylinder  where 


x 


(131) 


where  xa  is  the  axial  coordinate  and  LSHAPE  = 3 chooses  the 
values  for  a cone,  where 


x 


(132) 


These  formulas  have  been  found  to  be  quite  adequate  for  calcu- 
lating flows  over  long  bodies.  Small  initial  errors  in  the 
calculated  boundary  layer  become  negligible  in  a few  boundary- 
layer  thicknesses . 
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11.0  PROGRAM  OUTPUT 


11 . 1 STANDARD  OUTPUT 

Several  options  are  available  for  output  from  the  programs. 
The  parameter  NPRINT  on  item  number  2 in  Table  II  chooses  either 
all  the  available  output  (NPRINT  = 1)  or  only  that  essential  for 
monitoring  the  progress  of  an  iterative  calculation.  The  param- 
eter DXP  on  item  17  controls  the  printing  of  boundary- layer 
velocity  profiles.  For  iterative  calculations,  the  full  output 
from  the  inviscid-f low  program  is  printed  only  on  the  first 
iteration  (Nl  = 1).  For  subsequent  iterations,  only  the  short 
form  output  is  printed.  No  units  are  printed  in  the  output  since 
any  units  are  allowed  in  the  input. 

The  complete  program  output  is  presented  in  the  following 
list.  Steps  1-10  are  always  printed  for  the  first  iteration. 

Of  the  remaining  steps,  those  denoted  by  an  asterisk  (*)  are 
those  printed  at  the  end  of  an  xs-Bs  iteration  (N2  index)  when 
NPRINT  = 0.  For  NPRINT  = 1 all  output  (steps  11-18)  is  printed 
for  every  value  of  Nl. 

1.  Three-line  title  or  description. 

2.  List  of  all  values  of  integers  on  first  and  second  data 

card. 

3.  List  of  Body  Geometry  input. 

4.  List  of  other  input  values  for  inviscid  flow. 

5.  List  of  input  indices  for  boundary- layer  calculation. 

6.  List  of  Body  shape  data  for  boundary  layer. 

7.  List  of  other  boundary- layer  input  quantities. 

8.  Computed  geometric  parameters  in  normal  direction  for 

inviscid  flow. 

J - Normal  grid  index. 

AN  - Normal  coordinate. 

G - Stretching  function  derivative  (refs.  3 and  5). 

GH  - Stretching  function  derivative  at  half  intervals. 

9.  Computed  geometric  parameters  in  tangential  direction. 

I - Tangential  grid  index. 

S - Arc  length  along  reference  surface. 

X - Axial  coordinate. 

Y - Radial  coordinate. 
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THET  - Angle  of  reference  coordinate  surface,  6.  For 
closed  bodies,  0 is  the  same  as  the  body  angle, 
0B • For  open  bodies,  0 = 0b  on  the  forebody 
and  6 = 0 on  the  afterbody. 

Body  angle,  0g. 

Surface  curvature  on  closed  bodies.  For  open 
bodies  AK  is  the  surface  curvature  on  the  fore- 
body and  AK  = -(d2rw/dx^)  on  the  afterbody. 

Derivative  of  the  tangential  stretch  function 
(refs.  3 and  5) . 

10.  Plume  velocity  and  shape  distributions.  These  lists 
are  printed  at  the  beginning,  after  the  4th  iteration 
and  after  the  last  iteration  if  the  solution  converges 
to  the  least  squares  error  tolerance. 

11.  Inviscid  relaxation  iteration  history. 

IT  - Iteration  number. 

IT 

DPMAX  - Maximum  cb  correction,  max  d> . . 

il  J 

ID, JD  - I and  J location  of  DPMAX. 

RMAX  - Maximum  residual,  max  R. . I , where  R . . is 

ij  ill  il 

the  right-hand  side  of  the  difference 
equation. 

IR,JR  - I and  J location  of  RMAX. 

ISUB,ISUP  - Indicates  if  maximum  residual  occurred  at 
a subsonic  or  supersonic  point. 

RAVG  - Average  value  of  the  residual. 

RFl  - Relaxation  factor  for  subsonic  points. 

QF3  - Relaxation  factor  for  supersonic  points. 

NS  - Number  of  supersonic  points. 

SEC/CY  - Time  for  iteration  cycle. 

12.  List  of  solution  of  perturbation  potential. 

*13.  Tabulated  values  of  surface  pressure  coefficient,  Cp, 
Mach  number,  and  axial  velocity  on  the  body  along  with 
a rough  plot  of  C along  the  body.  This  plot  is 

r 


.IT-1 


THETB  - 
AK 

F 
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distorted  in  the  axial  direction  because  it  is  for  equal 
spacing  in  the  computational  space.  The  asterisks  show 

the  level  of  sonic  C . 

P 

*14.  Drag  coefficients  by  trapezoidal  integration  of  the  Cp's 
on  the  real  body.  The  displacement  surface  is  removed 
for  calculation  of  the  drag. 

*15.  Coordinates  x and  y of  the  sonic  line. 

16.  Boundary- layer  reference  velocity,  ue  , unit  Reynolds 

number,  R0  /L  and  viscosity,  . This  is  only  printed 

o 

on  the  first  step  (Nl  = 1) . 

*17.  List  of  boundary- layer  quantities  with  profiles  at 

intervals  governed  by  DXP.  This  is  always  printed  on 
the  4th  iteration  since  that  is  the  last  iteration  for 
which  the  calculation  begins  at  the  nose  of  the  body. 

AX  - Axial  distance  from  the  nose,  x. 

UTAU  - Friction  velocity,  uT . 

DELTA  - Boundary- layer  thickness,  6. 

DELST  - Displacement  thickness,  6*. 

i 

THETA  - Momentum  thickness,  0. 

CF  - Skin-friction  coefficient,  C^. 

UE/UZ  - Boundary-layer- edge  velocity  ratio;  ue/ue^. 

DELST+R  - Augmented  body  radius,  6*  + 

DX  - Integration  step  size. 

HTR  - Transformed  shape  factor,  6?/0^. 

*18.  Quantities  showing  status  of  inner  iteration. 

XMAX  - Location  of  maximum  change  in  boundary  layer 
6*  from  previous  iteration. 

DPMAX  - Maximum  change  in  6*. 

RBT  - Body  radius  at  nozzle  exit. 
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*19.  Rough  plot  and  list  of  u and  uc  . This  is  printed  at 

^ el 

the  end  of  the  inner  iteration  cycle  whenever  a new 
comparison  is  being  made  between  the  viscous  and  inviscid 
velocities.  This  is  followed  by  a list  of  the  current 
results  of  the  iteration. 

XSEP  = x 

s 

THET  = 0 (degrees) 

DRMS  = rras  error,  s (percent  of  u~  ) 

o 

11.2  SPECIAL  OUTPUT  MESSAGES 

Several  special  messages  are  contained  in  the  output  to  call 
attention  to  specific  conditions  that  may  occur.  The  messages 
are  listed  in  this  section  with  instructions  about  what  to  do 
when  they  are  encountered. 

(1)  DIVERGENCE.  RMAX  EXCEEDS  RCHEK, 

This  message  is  printed  by  the  inviscid-f low  program  if  the 
relaxation  procedure  diverges.  Check  all  input  to  verify  that 
it  is  correct.  If  no  obvious  errors  appear,  the  difficulty  is 
probably  either  in  the  choice  of  parameters  for  the  computational 
mesh,  or  the  smoothness  of  the  data  defining  the  body  shape. 

(2)  RF1  DECREASED  TO  BECAUSE  10-CYCLE  AVG  FOR  RMAX  INCREASED. 


This  message  refers  to  the  subsonic  relaxation  factor  in 
the  inviscid-f low  program.  The  initial  value  is  1.4.  The  value 
is  automatically  reduced  by  10  percent  if:  (1)  the  maximum 
residual,  averaged  over  10  cycles,  is  greater  than  that  for  the 
previous  10  cycles  and  (2)  the  last  maximum  residual  occurred  at 
a subsonic  point. 

(3)  QF3  INCREASED  TO  BECAUSE  10-CYCLE  AVE  FOR  RMAX  INCREASED. 


This  message  refers  to  the  supersonic  damping  factor  in  the 
inviscid-f low  program.  The  initial  value  is  0.1.  The  value  is 
automatically  increased  if:  (1)  the  maximum  correction,  averaged 
over  10  cycles,  is  greater  than  that  for  the  previous  10  cycles, 
and  (2)  the  last  maximum  residual  is  at  a supersonic  point. 

(4)  INPUT  FROM  TAPE13  HAS  INCOMPATIBLE  DIMENSIONS 

This  message  is  printed  if  the  dimensions  of  the  4> h j solu- 
tion read  from  Logical  Unit  13  (data  file  2 in  figure  3)  are  not 
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the  same  as  the  values  of  IMAX  and  JMAX  put  on  item  number  7 in 
Table  II. 

(5)  ****ITERATI0N  FOR  BOUNDARY  LAYER/INVISCID  FLOW  EQUILIBRIUM 
CONVERGED 

This  message  is  printed  whenever  the  maximum  change  in  6* 
between  iterations  is  less  than  the  specified  percent. 

(6)  METHOD  FOR  CALCULATING  UTAU  IN  DERIV  DOES  NOT  CONVERGE 

This  message  refers  to  the  iteration  used  to  solve  equation 
(55)  for  UT  when  6*  is  prescribed  in  the  boundary- layer  calcula- 
tion. The  only  known  cause  of  the  iteration  failing  to  converge 
is  an  error  in  the  input  data. 

(7)  DELTAI  HAS  BECOME  NEGATIVE 

STOP  INTEGRATION,  PRINT  PROFILE  AT  PREVIOUS  STEP 

This  message  refers  to  the  transformed  boundary- layer 
thickness,  6^.  The  error  condition  may  occur  due  to  the  initial 
integration  step  size  DXZ  being  too  large.  Another  possible 
cause  might  be  a too  sudden  change  in  the  body  shape,  or  in  the 
prescribed  ug  or  6*  distribution. 

(8)  DELTA  HAS  BECOME  NEGATIVE 

STOP  INTEGRATION,  PRINT  VALUES  AT  PREVIOUS  STEP 

This  message  is  not  expected  to  occur  in  the  finished  pro- 
gram. If  it  does,  check  the  input  data  carefully. 

(9)  METHOD  FOR  CALCULATING  K IN  ENTRAN  DOES  NOT  CONVERGE 

This  message  refers  to  the  iteration  used  to  solve  equation 
(72)  for  K when  6*  is  prescribed  in  the  entrainment- layer  calcu- 
lation. The  only  known  cause  of  the  iteration  failing  to  con- 
verge is  an  error  in  the  input  data. 

(10)  METHOD  FOR  CALCULATING  INITIAL  VALUE  OF  DELI  DOES  NOT 
CONVERGE 

When  initial  values  of  Cf,  6,  or  6*  are  known,  the  calcu- 
lation must  solve  an  integral  equation  for  the  initial  value  of 
the  transformed  thickness,  6j_.  This  is  done  by  iteration  in  a 
similar  manner  as  for  uT  described  in  message  (6) . If  the 
iteration  does  not  converge,  it  is  usually  due  to  errors  in  the 
input  quantities. 
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(11)  INTERMEDIATE  RESULTS  OF  ITERATION 

This  message  is  printed  at  the  end  of  each  step  of  the  xs-0s 
cycle.  It  is  followed  by  the  current  value  of  xs,  XSEP,  the 
current  value  of  es,  THET,  and  the  value  of  the  rms  error,  DRMS. 

(12)  ***FINAL  RESULTS*** 

BEST  SOLUTION  WAS  ITERATION  NO. 


This  message  is  printed  whenever  the  least  squared  error 
has  been  found  for  the  xs-0s  cycle  discussed  in  section  2.2. 

(13)  ****SOLUTION  CONVERGED  TO  LEAST  SQUARED  ERROR  TOLERANCE 

a 

This  message  is  printed  whenever  the  rms  error  is  less  than 

1 percent  of  uG  at  any  stage  of  the  iteration.  The  inviscid 
o 

and  viscous  solutions  immediately  preceding  this  message  are  then 
the  best  solutions  of  the  calculation  procedure. 

(14)  SKIN  FRICTION  HAS  BECOME  NEGATIVE  IN  AN  INCORRECT  MANNER. 
CHECK  ALL  INPUT  CAREFULLY 

This  message  will  be  printed  if  the  skin-friction  coefficient 
changes  sign.  It  may  indicate  that  the  initial  estimate  of  the 
separation  point  location  was  too  far  downstream.  It  has  usually 
been  observed  to  occur  when  strong  shocks  are  present,  or  when  too 
few  relaxation  steps,  MIT  were  used  initially. 

(15)  ATTACHED  BOUNDARY  LAYER  SOLUTION  TERMINATED  AFTER  

ITERATIONS 

This  message  will  be  printed  if  the  viscous-inviscid  inter- 
action converges  (Message  5)  and  no  separation  occurred  or  the 
separation  point  for  the  next  iteration  would  have  been  down- 
stream of  the  end  of  the  body. 


12.0  PROGRAM  OPERATING  PROCEDURE 

In  this  section,  the  construction  of  card  decks  for  opera- 
tion of  the  computer  programs  is  described.  First,  a general 
description  of  the  operations  required  is  given.  Then  the 
specific  Job  Control  cards  needed  for  operation  on  an  IBM  370 
computer  are  listed.  The  same  card  decks  should  be  applicable 
at  any  370  installation  with  minor  modifications. 
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12.1  GENERAL  JOB  CONTROL  SEQUENCE 

The  following  list  is  the  general  Job  Control  procedure  that 
would  be  required  to  run  the  programs  for  a complete  viscid- 
inviscid  interaction  calculation.  The  reader  is  referred  to 
figure  16  and  Table  I. 

1.  Create  partitioned  data  sets  for  restart  files 
(files  1-5  in  figure  16) . 

2.  Define  units  2 ,3, 8, 9,  and  10.  These  unit  numbers  are 
needed  for  output. 

3.  Define  units  11,12,13,14,  and  15  if  NRSTRT  = 1 in  the 
input  data.  These  unit  numbers  correspond  to  the  input 
files.  They  contain  data  created  in  a previous  run. 

For  starting  an  initial  calculation,  the  partitioned  data 
sets  would  be  created  in  a separate  operation.  Then,  since  no 
data  would  be  on  file,  only  units  2, 3, 8, 9,  and  10  need  to  be 
defined.  For  restarting  an  iterative  calculation,  all  data  files 
would  exist,  so  units  11  to  15  must  also  be  defined. 

To  execute  the  boundary- layer  program  alone,  unit  10  must 
be  defined  in  order  to  output  the  6 * + rw  and  ue/ueo  list.  Unit 
12  must  be  defined  when  LVARl  = 2,  and  unit  15  must  be  defined 
when  LDSTAR  = 1. 

To  execute  the  inviscid  program  alone,  units  2 and  8 must  be 
defined.  Unit  13  is  also  required  when  LREADP  = 1,  and  unit  14 
is  required  when  IXY  = 0. 

12.2  JOB  CONTROL  EXAMPLES 

In  this  section,  specific  examples  of  Job  Control  cards  used 
for  the  operations  discussed  previously  are  presented.  In  the 
examples,  the  computer  program  is  referred  to  as  "ITER"  with  the 
source  code  names  "SITER"  and  the  load  module  or  binary  version 
named  "BITER".  The  account  ID  used  in  the  examples  is  WYL.XM.K01. 
Logical  units  5 and  6 are  the  standard  input/output  file  numbers. 
It  is  not  necessary  to  specifically  define  these  unit  numbers  in 
the  JCL  deck. 

12.2.1  Creating  Partitioned  Data  Sets 

Partitioned  data  sets  for  use  as  input/output  disk  files 
must  be  created  before  the  normal  program  operation  can  proceed. 
The  following  procedure  is  suggested: 

Use  IBM  Utility  Program  IEFBR14 . 

Use  default  values  for  DCB  (DSORG=PO,RECFM=VS) . 
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On  3330  disk,  use  SPACE  in  tracks  as  follows  (refer  to 
figure  16  and  Table  I for  explanation  of  file  numbers) : 


VELBOD 

(File 

3) 

SPACE  = 

(TRK, (2,1,10) ) 

RESTRT 

(File 

5) 

SPACE  = 

(TRK , (10,2,10) ) 

PHI 

(File 

2) 

SPACE  = 

(TRK, (20,4,10) ) 

XOFILE 

(File 

1) 

SPACE  = 

(TRK,  (4,1,10) ) 

DSFILE 

(File 

4) 

SPACE  = 

(TRK,  (6,1,10)  ) 

Example  of  creating  a partitioned  data  set  called  VELBOD: 
//EXEC  PGM=IEFBR14 

//A  DD  DSN=WYL. XM. KOI . VELBOD ,V0L= volume , 

//  UNIT=333Q, DISP= ( ,CATLG)  , 

//  SPACE=(TRK, (2,1,10) ) 


12.2.2  Starting  an  Iteration  Sequence 

To  start  an  iteration  sequence,  unit  numbers  2,  3,  8,  9, 
and  10  must  be  defined.  The  specific  sequence  of  cards  used  to 
perform  the  calculations  presented  in  section  13.1  is  as  follows: 


//EXEC  FORTGO , P ROG= ITER , VOL=VO lume 
//  LIB='WYL.XM. KOI. BITER' 

//GO. FT02F001  DD  DSN=WYL. XM. KOI .VELBOD (RUNl ) , 
//  DISP=OLD 

//GO. FT03F001  DD  DSK=WYL . XM. KOI . RESTRT (RUNl) , 
//  DISP=OLD 

//GO.FT08F001  DD  DSN=WYL. XM. KOI . PHI (RUNl) , 

//  . DISP=OLD 

//GO.FT09F001  DD  DSN=WYL. XM. KOI .XOFILE (RUNl) , 
//  DISP=OLD 

//GO. FT10F001 ' DD  DSN=WYL. XM. KOI , DSFILE (RUNl ) , 
//  DISP=OLD 

//GO. SYS IN  DD  * 


( 


Input  data  cards 


/* 
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12.2.3  Restarting  an  Iteration  Sequence 

The  specific  cards  used  to  perform  a restart  of  the  calcu- 
lation started  in  the  previous  section  are: 

//  EXEC  FORTGO/PROG=ITER,VOL= volume, 

//  LIB='WYL.XM. KOI. BITER' 

//GO . FT 0 2F 0 01  DD  DSN=WYL. XM. KOI .VELBOD (RUN2)  , 

//  DISP=OLD 

//GO. FT03F001  DD  DSN=WYL. XM. KOI . RESTRT (RUN2) , 

t 

//  DISP=OLD 

//GO . FT 0 8F 0 0 1 DD  DSN-WYL . XM.K01 . PHI (RUM 2)  , 

//  DISP=OLD 

//GO.FT09F001  DD  DSN=WYL. XM. KOI. XOFILE (RUN2) , 

//  DISP=OLD 

//GO.FT10F001  DD  DSN=WYL. XM. KOI .DSFILE (RUN2 ) , 

//  DISP=OLD 

//GO . FTllF 001  DD  DSN=WYL. XM. KOI. RESTRT (RUNl) , 

//  DISP=OLD , LABEL= { , , , IN ) 

//GO. FT12F001  DD  DSN=WYL. XM. KOI .VELBOD (RUNl) , 

//  DISP=OLD,LABEL=(, , ,IN) 

//GO . FTl 3F 0 01  DD  DSN=WYL. XM. KOI .PHI (RUNl ) , 

//  DISP=OLD,LABEL= ( , , , IN) 

//GO.FT14F001  DD  DSN-WYL. XM. KOI .XOFILE (RUNl) , 

//  DISP=OLD , LABEL= ( , , , IN ) 

//GO.FT15F001  DD  DSN=WYL . XM . KO'l . DSFILE  (RUNl)  , 

//  DISP=0LD , LABEL= ( , , , IN } 

//GO. SYS IN  DD  * 


C 


Input  data  cards 


r 


These  cards  were  used  with  the  example  discussed  in  section  13.2. 
The  input  data  cards  required  for  restarting  are  summarized  in 
Table  III. 
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12.2.4  Executing  the  Boundary-Layer  Program  Alone 

The  specific  cards  used  to  perform  the  calculations  dis- 
cussed in  section  13.3  are  listed  in  this  section.  In  the  exam- 
ple shown  here,  all  input  is  assumed  to  be  from  cards,  but  the 
output  list  of  6*  + r„  and  ue/ueQ  is  to  be  saved  on  unit  10. 

Unit  12  would  be  required  for  input  if  LVAR1  = 2,  and  unit  15 
would  be  required  if  LDSTAR  = 1.  The  cards  used  in  the  example 
in  section  13.3  are: 

//  EXEC  FORTGO , PROG=ITER,VOL=volume 
//  LIB=TWYL.XM. KOI. BITER' 

/ / GO . FTl OF 0 0 1 DD  DSN=WYL. XM. KOI . DSFILE (RUNl)  , 

//  DISP=0LD 

//GO. SYS IN  DD  * 

^ Input  data  cards 

/* 


12.2.5  Executing  the  Inviscid  Program  Alone 

The  cards  used  to  perform  the  calculations  discussed  in 
section  13.4  are  listed  in  this  section.  In  this  example,  the 
velocity  potential,  c|>,  is  input  from  unit  13,  and  the  new  solu- 
tion for  $ is  output  on  unit  8.  The  calculated  velocity  on  the 
body  is  output  on  unit  2.  Input: from  unit  13  corresponds  to 
LREADP  = 1 in  the  card  input  data.  In  addition,  unit  14  would 
be  required  for  input  of  the  body  shape  if  IXY  = 0 in  the  card 
input  data.  The  specific  cards  used  in  the  example  are: 

//  EXEC  FORTGO, PROG=ITER,VOL= volume, 

//  LIB='WYL.XM. KOI. BITER' 

/ / GO . FT 0 2F 0 0 1 DD  DSN=WYL. XM.K01 . VELBOD (RUN2 ) , 

//  DISP=OLD 

//GO.FT08F001  DD  DSN=WYL. XM. KOI .PHI (RUN2) , 

//  DISP=QLD 

//GO .FTl 3F 001  DD  DSN =WYL. XM. KOI . PHI (RUNl) , 

//  DISP=OLD , LABEL= ( , , , IN ) 

//GO. SYS IN  DD  * 

^ Input  data  cards 

/* 
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13.0  NUMERICAL  EXAMPLES 

In  this  section,  several  example  calculations  are  presented 
to  aid  in  program  checkout.  An  example  is  presented  of  a com- 
plete viscid-inviscid  interaction.  An  example  is  also  presented 
of  the  use  of  the  boundary- layer  program  alone  for  a two-dimen- 
sional geometry.  That  example  also  demonstrates  the  two  options 
for  boundary  conditions,  having  ue  specified  in  the  beginning  of 
the  calculation,  and  6*  specified  in  the  second  part.  Input  data 
for  another  sample  case  are  also  presented  to  demonstrate  the 
use  of  the  program  to  calculate  the  inviscid  flow  alone. 

13.1  AXISYMMETRIC  INTERACTION 

A list  of  the  punched  card  input  data  for  a sample  calcula- 
tion on  the  boattailed  body  shown  in  figure  8a  and  b is  presented 
in  figure  17.  The  case  being  calculated  is  for  a free-stream 
Mach  number  of  0.8.  The  body  corresponds  to  the  ogive-cylinder 
body  with  a circular-arc  boattail  described  in  reference  20.  The 
JCL  card  deck  for  this  case  has  been  presented  in  section  12.2.2. 
The  running  time  for  the  complete  calculation  is  about  3.6  minutes 
on  the  IBM  370/165. 

Selected  output  for  the  sample  case  is  shown  in  figure  18. 

The  solution  converged  to  the  least  squared  error  tolerance  in 
25  iterations. 

The  complete  list  of  output  for  this  case,  with  NPRINT  = 0, 
consisted  of  a total  of  approximately  2500  lines.  Output  steps  - 
1-7  as  listed  in  section  11.1  have  been  omitted  from  this  presen- 
tation since  they  simply  verify  the  input  data.  The  output  pages 
shown  are  those  corresponding  to  steps  8-19  of  the  set  described 
in  section  11.1  for  NPRINT  - 0 'on  the  first  iteration  and  then  a 
few  samples  of  the  rough  plots  of  u^  and  ue-j-  at  intermediate 
steps,  including  that  at  iteration  25  {fig.  18(f)).  Finally, 
figure  IB  concludes  with  the  inviscid  solution  and  boundary- layer 
and  plume  solutions  corresponding  to  the  final  result. 

13.2  EXAMPLE  OF  RESTARTING  AN  INTERACTION  CALCULATION 

The  punched  card  input  data  for  restarting  the  interaction 
calculation  of  section  13.1  if  that  calculation  should  undergo 
premature  termination  at  an  intermediate  iteration  is  presented 
in  figure  19.  The  JCL  card  deck  for  this  calculation  was  pre- 
sented in  section  12.2.3.  Note  that  the  iterative  calculations 
cannot  be  restarted  at  any  arbitrary  iteration  using  the  restart 
file,  unit  11.  That  file  and  the  other  output  files  contain 
only  the  data  that  were  output  just  prior  to  the  termination. 
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13.3  TWO-DIMENSIONAL  BOUNDARY  LAYER 

A list  of  the  punched  card  input  data  for  a sample  calcula- 
tion on  the  two-dimensional  configuration  shown  in  figure  20  is 
presented  in  figure  21.  Note  that  two  sets  of  input  are  pre- 
sented for  this  case,  giving  an  example  of  the  options  of  pre- 
scribed u„  and  prescribed  6*.  The  data  for  ue  and  6*  were 
obtained  from  the  experimental  results  of  reference  18  which 
indicate  separation  occurring  in  an  adverse  pressure  gradient 
region  downstream  of  a shock  wave.  The  output  for  the  complete 
boundary-layer  calculation  are  presented  in  figure  22.  No 
external  data  files  were  used  for  input  for  this  case.  The  JCL 
card  deck  for  this  case  was  presented  in  section  12.2.4. 

13.4  AXI SYMMETRIC  INVISCID  FLOW 

The  punched  card  input  data  for  a sample  calculation  of  the 
inviscid  flow  alone  are  presented  in  figure  23.  The  output  for 
this  case  is  shown  in  figure  24. 


i 
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Figure  2.-  Schematic  af  viscous- inviacid  interaction  iteration. 


74 


Figure  4.-  Boundary  layer  coordinates  and  notation. 
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Figure  4.-  Boundary  layer  coordinates  and  notation. 
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S P 

(c)  Friction  velocity  and  displacement  thickness. 


Figure  5.-  Typical  variation  of  boundary  layer 
quantities  for  separated  boattail  flow. 
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Figure  6.-  Schematic  of  exhaust 
flow  entrainment  me 
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plume- separated 
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Figure  7.-  Typical  velocity  profiles  in  the 
entrainment  region. 
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(a)  Forebody, 


All  dimensions 
are  referenced 
to  the  maximum 
diameter 


(c)  Boattail  configuration  2. 

Figure  8.-  Body  with  circular  arc  boattail  (ref.  20). 


(b)  M =0,8. 
o 


Figure  9.-  Comparison  between  calculated  and  measured 
pressure  coefficient  distributions  on  unseparated 
boattail- sting  configuration. 
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Figure  10.-  Comparison  between  calculated  and  measured 
pressure  coefficient  distributions  on  a boattail- sting 
configuration  with  separation. 
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(c)  Mo  «*  0.9. 


Figure  10.-  Concluded 


(c)  Skin-friction  coefficient 


Figure  11.-  Concluded. 
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X/D 


(c)  Ptj/P0  = 2,0>  Mo  = °-6* 


Figure  13.-  Concluded. 


1.0 


2.0 


4.0 


3.  0 


/P 


o 


(b)  M = 0.6. 
o 


- Comparison  of  experimental  and 
.ail  drag  on  configuration  1 with 
simulated  with  high  pressure  air 
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Figure  16.  General  relationship  of  programs  and  data  files. 
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NASA  DATA  COMPARISON  - CONFIGURATION  1 
JET  SIMULATED  WITH  HIGH  PRESSURE  AIR 
BOATTAIL  L/D  = 0.8  DE/t)  =0.51 

000  700  0 

0 11110 


1.4 

98 

0.  8 

0.0 

0.0 

0. 1 

0.025013 

0.  2 

0.050026 

0.3 

0.075038 
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0.100051 

0.5 
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0.6 

0. 150077 
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0. 225115 
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1.  1 
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0.300154 
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1.4 
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1.  5 

0. 375192 

1.6 

0.400205 

1.  7 

0.425218 

1.  8 

0.450230 

1.9 

0.473630 

2.0 

0. 489760 

2.  1 

0.498366 

2.  2 

0.500000 

2.3 

0. 500000 

2.5 

0.500000 

3.0 

0.500000 

3.5 

0. 500000 

4.0 

0.500000 

5.0 

0.500000 

6.0 

0. 500000 

7.0 

0.500000 

7.5 

0.500000 

8,0 

0. 500000 

8,025 

0.499781 

8.05 

0.499125 

8.075 

0.498030 

8.  1 

0.496496 

8. 125 

0.494521 

8.  15 

0.492104 

8. 175 

0.489241 

8.  2 

0.485931 

8.  225 

0.482171 

(a)  First  49  cards. 

Figure  17.-  Input  data  for  viscous-inviscid 
interaction  calculation. 
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8.  25 

0.477956 

8.  275 

0.473282 

8.3 

0.468146 

8.325 

0.462542 

8.35 

0.456463 
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8.6 

0. 367897 

8.625 
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8.65 
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0. 302368 

8.  75 
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8.775 

0.271518 

8.  8 

0. 255000 

8.  81 

0. 255000 

8.  82 

0. 255000 

8.  83 

0. 255000 

8.  84 

0. 255000 

8.  85 

0. 255000 

8.  86 

0. 255000 

8.87 

0. 255000 

8.88 

0. 255000 

8.  89 

0. 255000 

8.  9 

0. 255000 

8.92 

0.255000 

8.94 

0. 255000 

8.  96 

0. 255000 

8.98 

0. 255000 

9.0 

0. 255000 

9.025 

0. 255000 

9.05 

0.255000 

9.075 

0. 255000 

9.1 

0. 255000 

9.  15 

0. 255000 

9.  2 

0. 255000 

9.  3 

0. 255000 

9.4 

0. 255000 

9.5 

0. 255000 

9.7 

0. 255000 

9.9 

0. 255000 

(b)  Next  49  cards. 
Figure  17.-  Continued. 
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10.1  0.255000 

10.3  0.255000 

10.5  0.255000 

10.8  .0. 255000 

11.2  0.255000 

11.6  0.255000 

12.0  0.255000 

61  31  0 0 0 
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0.0 
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» 

o 
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o 

8.8 

1.0 

8.  0 

T 
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1.4 

4.0 
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0.0 

523.8 
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0.0125 

F 

1 

12  3 

1 0 

0 

0.  5 

2134. 1 

576.9 

0.97 

0.  15 

10,0 

0. 166667 

0.0 

(c)  Remaining  15  cards. 
Figure  17.-  Concluded. 


95 


AEDC-TR-79-4 


nopm\l  coord.  *;thftch  for  ai_f=  1.300 


J AN  <;  GH 


J 

0.349CF- 

-5-3 

0.0 

0 .7T25E-94 

2 

0 ,16171- 

. P3. 

0_.  154  55.-03 

0.4fc07F-03_  

3 

0 • 6 3395 

0? 

0-7C6ar-03 

0.1  3£  bF-02 

4 

0 . 360 at- 

02 

0 . 196  4E— 02 

0 .2900E-C2 

5 

0*2  39  IF 

02 

0.3836F-02 

0.51486-02 

. . 6 

0 .17206 

0.2 

C.6460E-02 . 

0.81 82.E  — 02 

7 

0,*  3 0?F 

02 

0.990SC-02 

0 • 1 207E-r 1 

0. 

.QiJL-Q  2 2f. 

0,  l*24fc-Q \ 

<Ui  687^0A_  — . 

9 

0-821  5h 

01 

0 -L951E-01 

0.2  2f  5C-01 

10 

0.*729E 

01 

0.?5795~01 

O.2947E-.01 

1 1 

0.5S88C 

0 1 

0.331  45-01 

O.J730E-OI 

.12 

0-469  06 

01 

0.4  16  1C— 0 1 

0.46446-01 

1 3 

0 • "*96  QE 

01 

9*SJ  ?6F1— 0 1 

0.5670E-01 

14  . 

U«13Z7r_2_i. 

_QjJLL1  &iir  QJL  _ 

.0-6  05.4^T0J 

1 5 

o.23ae»i 

n * 

0.7  42  3r-0 1 

0.93 1 OF  — 0 1 

- - 4L 

0.2.474F 

OJ 

.&.F7^7F-01  . 

0.9  535E-01.  _ .. 

17 

9.2123C 

01 

C.1Q28C  00 

0-1  1 1 OF  00 

_ . 18 

.0 . 1 322F 

.0.1 

0..  1193C  00 

0.1282F  0.0... 

19 

0.!  5C  IF 

01 

0.13721*  00 

0.1  4 7 OF  00 

rif.13  3.aC!. 

o.o  . 

0 *X_6_74E  _ OQ 

2 1 

0. 1 139E 

01 

0.1 78  OE  OO 

0.1895F  00 

22 

0-9585F_ 

_00 

_ 0 .2  01  OF  00 

0.21  J4t  _0O 

23 

0 .80  20= 

00 

0.2259E5  00 

0.P391E  00 

2 4 

0 .662 ?fc 

no 

n.25£4E  00 

0.2  66  76  00. 

25 

0 .537 ie 

00 

0.28105  00 
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(a)  Computed  geometric  parameters  in 
normal  direction  for  inviscid  flow. 

Figure  18.-  Selected  output  for  viscid-inviscid  integration  calculation. 
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Figure  18.-  Continued. 
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Figure  18.-  Continued 
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Figure  IB.-  Continued 
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..  i_,432er 

03 

_ . _ 1 #210  4T 

0 1 

2.9245E-Q 1 

1 ,213  4£ 

01 

1 * 4 3285 

03 

1 ,23ft4£ 

01 

"~2  ,9 245E— 0 1 

1 ,238*C 

01 

1 ♦432BE 

03 

1 .294 4E 

01 

. 2 • 9245E -0 1 

1 .2944E 

01 

1 .432  MI 

03 

1 .35Q4E 

01 

2,O245E-01 

1 ,3  504K 

01 

1 .4320F 

03 

.1,40645 

01 

2.9246jF-pl 

1 « 4 Ofi  4E 

01 

1 ,43285 

03 

1 ,462  45 

o i 

" ""*  2,924  SE - 0 1 

l ,4f 24E 

01 

1.432SE 

03 

1 ,630  4E 

11 

2 ,924  55 -0^ 

1 •63A4EI 

01 

1 .432HE 

03 

1 ,7904E 

91 

2,9245E-oT 

1 ,798  45 

01 

] -432BF 

o? 

1.9663F 

01 

2,9$4SE-01. 

o 

SO 


(i)  Concluded. 
Figure  IS.-  Concluded. 
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AEDC-TR-79-4 


NASA  CONFIGURATION  1 
(BLANK  CARD) 

RESTARTING  ITERATION  SCHEME  AFTER  ITERATIONS  FOR  RUN 
1 0 0 600  0 


Figure  19.-  Input  data  for  restart. 
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ALBER  TEST  CASE 

TWO-DIMENSIONAL  BOUNDARY  LAYER  CALCULATION 
SPECIFYING  VELOCITY  DISTRIBUTION  FROM  X=0 


0 

1 

1.4 

1 

22 


0 

0 


0 0 

0 0 

. 7325 
0 0 


0 

2 


0 


0.  0 

9750. 

1.0 

9480. 

2.0 

9400. 

3.0 

8500. 

4.0 

8650. 

5.0 

9750. 

6.0 

11650. 

7.0 

12980. 

8.0 

14200. 

9.0 

15480. 

10.0 

16100. 

10.  25 

14800. 

10.  75 

13400. 

11.0 

12600. 

11.50 

12100, 

12.0 

11150. 

13.  0 

10900. 

14.0 

10800. 

15.0 

10300. 

16.0 

9900. 

17.0 

9800. 

18.0 

10200. 

1.0 

15.0  585.0 

0.  0 

15.0  -5.0 

.00233 

.065 

19 

0.0 

0.0 

1.0 

0.0 

2.0 

0.0 

3.0 

0.0 

4.0 

0.  25 

5.0 

0.50 

6.0 

0.70 

7.0 

0.85 

8.0 

0.98 

9.0 

1.0 

10.0 

0.98 

11.0 

0.  85 

12.0 

0.70 

13.0 

0.  50 

14.0 

0.  25 

15.0 

0.0 

16.0  ’ 

0.0 

17.0 

0.0 

18.0 

0.0 

(a)  Input  for 

Figure  21.-  Inp 

.99 

0.0 


0.0 

0.0 


0.  0 
0.0 


0.  0 
0.0 


0.0 


ue  specified. 
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ALBER  TEST  CASE 

TWO-DIMENSIONAL  BOUNDARY  LAYER  CALCULATION 
SPECIFYING  DELST  DISTRIBUTION  FROM  X=10. 75 

00000 

1 0 0 0 2 0 

1.4  .7325 

2 0 0 0 2 0 0 

13 

10.0  .03 

10.25  .05 

10.75  .07 

11.0  .083 

11.5  .112 

12.0  .148 

13.0  .314 

14.0  .435 

14.5  .460 

15.0  .440 

16.0  .306 

17.0  .230 

18. 0 . 190 

1.0  15.0  585.0  .99  0.0  0.0  0.0  0.0 

10.75  18.0  -5.0  0.0  0.0  0.0  0.0 

.00077 

13400.  -2820. 

19 

0.0  0.0 

1.0  0.0 

2.0  0.0 

3.0  0.0 

4.0  0.25 

5.0  0.50 

6.0  0.70 

7.0  0.85 

8.0  0.98 

9.0  1.0 

10.0  0.98 

11.0  0.85 

12.0  0.70 

13.0  0.50 

14.0  0.25 

15.0  0.0 

16.0  0.0 

17.0  0.0 

18.0  0.0 


(b)  Input  data  for  case  with  6*  specified. 


Figure  21.-  Concluded. 


1 13 


UZ  * 9.90272E  03 
R 6Z  - J.02425E  05 
NU  Z = 3.27A43E-02 


AX 

UTAU 

DELTA 

DELST 

THETA 

CF 

UE/UZ 

DELST *R 

OX 

HTR 

J«0 

3.32242C 

02 

3.901 40E-01 

6.52716E-02 

3.68389E-02 

2.  33000E-03 

0.98458 

6.52716E-02 

0.0 

1 .1853 

5,000  OOE- 

-01 

3.21310b 

02 

3 . 958  C3E- 0 1 

6.93925E  — 02 

3.912S3E-Q2 

2.2470SE-D3 

0.97097 

6.939Z5E-02 

1.250E-01 

1 .2071 

1 «ooooot 

00 

3.1053 5E 

02 

4. 10471E-01 

7.37869E-02 

4.153 99E-02 

2.16527E-03 

0.95731 

7.37B69E-C2 

1.250E-01 

1.2288 

i .bojggr 

00 

3.0ei82E 

02 

4.201  19E-01 

1 . 5554  DC  —02 

4.258836-02 

2. 15244E-03 

0.95326 

7.S5540E-02 

1.250E-01 

1.2293 

2 « 000  00  £ 

00 

3.058ICE 

02 

4.298C8E-01 

7.  7JS37K-02 

4.36547E-02 

2.13935E-03 

0.94923 

7.73S37E-02 

1 .250E- 01 

1 .2300 

2-bOOOOC 

00 

2*6746 6E 

02 

4.51453E-01 

9.*03*2E-02 

5.1 9691 E -02 

1 . 82001 E— 03 

0 ,90409 

9.40  342E-02 

1.250E— 01 

1.3310 

3.00000C 

00 

2.2596?f 

02 

4.654gaE-01 

1.  1805  IE—  0 1 

6.2S7S7E-02 

I ,45 35 8 E-0  3 

0.85835 

1 . 1BOS1E-01 

1.250E-01 

1 .*640 

3.500  00b 

00 

2. 38  37  7 E 

02 

4 • 91889E- 0 1 

la  12308E— 01 

6. 086 59 E- 02 

1 .50704L-03 

0.86543 

2.37 J08E-01 

1 .250E-0L 

1.4072 

4 .ooo  oot: 

00 

2. 4909*E 

02 

4.96767E-01 

1 .076776-01 

5.93432E-02 

1.70071 E-03 

0.07350 

3.57677E-0I 

I.250E-01 

1.3613 

4.50000E 

00 

3.00  00  4E 

02 

4.8  79  66  E—  0 1 

a.404 iae-02 

4.61035C-02 

2. 15065b— 0 3 

0.92949 

4.590426-01 

1 .2SDE-01 

1 .2059 

5.00C OJE 

00 

3.43  19*£ 

02 

4.89869E-01 

6.543496-02 

3 ■ 9 39  73  E— 02 

2.48613E-03 

0.98458 

S.69435E-01 

1.250E-01 

1 .098* 

5.5  00  00E 

00 

4j».l  2 IB  4£ 

02  . 

S . 1 268 /£ — 0 l__ 

5.259962-02 

2.6856aE-i? 

2.9Q462E-Q3 

1_.0  8214_ 

6.S2603E-01 

1 .25OE-01 

0.9467 

6. ODD  ODE 

00 

4.72248E 

02 

5.79769E-01 

4.3223SE-02 

2.1 19696-Q2 

3.  15095E—03*" 

1 . 1 7644 

7. 432236-01 

1.250E-01 

0.8087 

<>•500  OOE 

00 

5.  1083  ir 

02 

6 *996 77b- 01 

3.S5237E-0 2 

1 .650496-02 

3. 23355E-03 

1.24456 

8. 14524E— 0 1 

1 .2505—  01 

0.6903 

7*000  00E 

00 

5. 46276E 

02 

4.733Q4E-01 

3.81 249E-02 

1 . 73663E-0? 

3. 26948E— 03 

1.3107S 

8.881256-01 

1.250E-01 

0.8184 

7.500 DOL 

00 

5.84  73  7b 

02 

5.411 B9E-01 

3.61652E-02 

1 .457276-02 

3. 34613E— 03 

1.37327 

9.511 6SE-01 

1.2S0E-01 

0.7332 

8 , 000  ODE 

00 

6. 15961b 

02 

7 .38021 E—  0 1 

3.61 176E-02 

1 .10344E-02 

3.33669E-03 

1.43395 

1.01612E  00 

1 .250E-OI 

0.59S4 

8.5  00  OOE 

00 

6.63 13  BE 

02 

3.24J 66E-01 

3.25074E-02 

I.35469E-02 

3.45291 E-03 

1.49974 

1.02251E  00 

1 .250E-01 

0.8748 

9.DOOOOC 

00 

7.08463L 

C 2 

3.5 1 y 5SE- Q 1 

3+195/7E-02 

1.22L38E-02 

3. 54016E-O3 

1 .56321 

1.031 96E  00 

1 *250E- 01 

0.821 l 

9,500 00b 

00 

7.2888  It 

02 

J.56377E-01 

3. 252266-02 

1 .20838  E— 02 

3.55447E-03 

1*59481 

1.02252E  00 

1 .250E-01 

0.8153 

i voooooc 

01 

7.49256L 

02 

3.609 C2E- 01 

3.3I672E-Q2 

1 *199006-02 

3.S6767E— 03 

1.62581 

1*0131 7E  00 

1.250E-0I 

0.8107 

1 .050  ODE 

01 

5.69590E 

C 2 

1.97371 E—  0 1 

3.5104SE— 02 

1.62790E-02 

2. 09 74 OE— 03 

1.42514 

9+50 1 04E— 01 

1 * 250E-  01 

1.1481 

1 *100  ODE 

01 

4.02644E 

02 

2. 013486-01 

S.  125616-02 

2 .364866—02 

1.90670E-03 

1.27230 

9.01256E-01 

1.Z50E-OI 

1 .4164 

i .ISO  OOE 

01 

3.53120E 

02 

2. 17B10E-01 

6.C16S4E-02 

2.76030 E— 02 

i .61324E-03 

1 .221 08 

8.35165E-01 

1 .2SOE-01 

t .5065 

1 .200  OOE 

01 

2. 34226E 

02 

2.b649DE-0l 

8.96149E-02 

3.671 89E-02 

8.5721  IE— 04 

1.12595 

7.896 158-01 

1.250E-01 

1 .8657 

1 .250  ODE 

01 

2.49992E 

02 

2*715 78 E~ 01 

8.52057E-02 

3.83132E-02 

1.00 1 S2E-03 

1*11336 

6*  89? OOE— 0 1 

1.230E-01 

1.7627 

1 .300  03c 

01 

2.57 45 2E 

02 

2. 8 7398E- 01 

9.0540ZE-OZ 

4*0  05506—02 

l. 09042E-03 

1.10071 

5.905406-01 

1.2S0E-0I 

1.7029 

1 ,J50  03 E 

Cl 

2.71775E 

02 

2.00251 E— 0 1 

8.  91  909E-02 

4.08723E-02 

l*  22  784003 

1*09566 

4.64191E— 01 

1 . 250E-  0 1 

1.62*3 

1.400 OOE 

01 

2.81  560b 

02 

3.13199E-01 

8.  £76286-02 

4.172 75E-02 

1.33169E-03 

1.09061 

3.38763E-01 

1.250E-01 

1 .5686 

1 *450  OOE 

01 

2.6290 3E 

02 

3.34069E-01 

9.736*46-02 

4 *543696-02 

1 .2237 1 E-03 

1.06547 

2.22564= -0  t 

1 .250E-01 

1.6110 

1 *50  3 OOE 

ex 

2.4358 5E 

02 

3.57543E-01 

1 . C7992E— 0 1 

4.96323E-02 

1 . 10873E-03 

1.04012 

1.D7992E-01 

1.250E-01 

1.6616 

(a)  Output  from  case  with  specified. 

Pigure  22.-  Output  from  two-dimensional  boundary- layer  calculation. 
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uZ  = 9.9027ZF  03 
rl£S  = J.0242I.L  OS 
NUZ  ■ J.?Z«4Jt-02 


AX 

UTAu 

DEL  TA 

DkLST 

Theta 

1 .U  75  006 

01 

2. r&asjE 

02 

i . raa  7<je-oi 

7. Ot  000k -02 

2 .497736—02 

1 

01 

2.  40  57 2 E 

02 

1 .366646-01 

7.64999E-02 

2 .643906-02 

1*130  0 0 J- 

01 

2 *2 1972k 

02 

1 .9662  U--0  1 

0. 299996-02 

2.795126-02 

1 . l ‘>0  JOc 

01 

1 . 46  28 OE 

02 

dm  4565BE-  01 

1 . 12030E-01 

3. 39O73E-02 

i .200001: 

01 

9.  lb3S3t 

01 

3.  046  55 E—  0 l 

1 . *80  00k  — 0 1 

4.039236-02 

i *2'jJ  3JE 

OL 

-y.2o©3  0E 

00 

4.2396  IE— 0 1 

2.31 OOOE  — 01 

4.79995C-02 

I.30JUJL 

01 

- S • 4 2 52  3E 

01 

S. 44/06E-01 

3. 140COE-0 t 

5.361576-02 

1 - J50 03  E 

01 

"(is  36  21  4£ 

ot 

6.38-)68k-01 

3 ♦ 74500k -01 

5.951156-02 

1 .4  00  00k 

01 

-7.71 28 2 £ 

0 1 

7.291  38001 

4.3S000E-01 

6 • 358 75 E— 02 

1 .43000c 

01 

-7m  ld5o4k 

01 

T • 76346E-01 

4.600006-0 1 

6- 9 2 785 E— 02 

I .50300k 

01 

-4.45844E 

0 1 

7 .735246*-  01 

4. 400  00k  — 0 L 

7.833586-02 

1 • 330  00c. 

01 

2.0022dE 

01 

7.1 S6Cje-01 

3 • 730C0E— 0 1 

8.69426k— 02 

1 .feOMOOl 

01 

6.27963k 

0 l 

6.36124E— J1 

3. C6000E-01 

9. 0 543  1 E- 02 

1 .0  50  00  C 

01 

1. 25859k 

02 

6.34677k- 01 

2. 630  00k- 0 1 

9.  1 53196-02 

1 .700  00  k 

31 

1. 7 269  4E 

02 

6.  1 26  1 7fc-  01 

2.300006-01 

8.942 ftfeE-02 

1 . 7ti0  00i: 

01 

2.0264  1C 

02 

6.06530k— 0 1 

2.  100006-01 

8. 760846-02 

1 .HOC  00 E 

01 

2. J417oC 

02 

6. 034 1 IE- 01 

1.9000 Ok -01 

8.4  36  94  E — 02 

CF 

ijc/ui 

OEUST+R 

OX 

HTR 

7.71 98  IE  — 04 

1.35316 

9.52500E-01 

0.0 

1.9901 

6.407006-04 

1.33210 

, 9.42750E-O1 

3.125E-02 

2.0924 

5. J7044E-0* 

1.31357 

9.33000k- 0 1 

3. 125E-02 

2.1881 

2. 60983E-04 

1.25296 

0 » 870006-0} 

1 .2506-01 

2.5641 

1 • I2095E-04 

1 .20558 

8*  48000k- 01 

1.250E-OI 

2.9561 

-9.75670C-O7 

1.17132 

8.31 OOOE-Ot 

l .250E-01 

4.0823 

—4.3331 5E-05 

1.15562 

8.140OOE-OI 

1 .2506-01 

5.1393 

-6.518816-05 

1.14165 

7. 495006-01 

1 .250k- 01 

5- 6144 

—9.1 3034k— 05 

1.13460 

6.850006-01 

1.250E-01 

6.2043 

-B. 1 1 766E-05 

1.12291 

5 * 8500  OE- 0 1 

1 .250k-  01 

6.0171 

-3. 262696-05 

1.10162 

4.40  OOOE-Ot 

1 .250k—  01 

4.9922 

6.884646-06 

1.06002 

3.73000E-01 

1 .250E-  01 

3.7Q74 

1.20372E-0* 

1.06939 

3. 06000k- 01 

1 .250k- 01 

2.8259 

2.80038E-04 

1.06617 

2.63000E-0! 

1 .250k- 01 

2.3087 

5*  2I237E-04 

1.07271 

2.30000E-01 

1 .250h-0l 

2.034? 

7.082576-04 

1.076U6 

2. ! JOOOE-Ol 

1.250E-01 

1.8524 

9. 19959E-04 

1.00125 

l .90000E-01 

t .2  506-01 

1 .6961 

tb) 


Output  from  case  with 


f*  specified. 


Figure  22 


Concluded, 
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NASA  DATA  COMPARISON  - CONFIGURATION  1 
JET  SIMULATED  WITH  A SOLID  STING 
BOATTAIL  L/D  -0.8 

000  700  0 


-1 

0 0 1 

1.4 

98 

0.9 

0.0 

0.0 

0.  1 

0.025013 

0.  2 

0.050026 

0.3 

0.075038 

0.4 

0.100051 

0.5 

0. 125064 

0.6 

0. 150077 

0.7 

0. 175090 

0.8 

0. 200102 

0.9 

0. 225115 

1.0 

0. 250128 

1.1 

0.275141 

1.2 

0. 300154 

1.3 

0.325166 

1.4 

0. 350179 

1.5 

0.375192 

1.6 

0.400205 

1.7 

0.425218 

1.  8 

0.450230 

1.9 

0.473630 

2.  0 

0.489760 

2.  1 

0.498366 

2.  2 

0. 500000 

2.  3 

0. 500000 

2.5 

0. 500000 

3.  0 

0. 500000 

3.  5 

0. 500000 

4.  0 

0. 500000 

5.0 

0. 500000 

6.0 

0. 500000 

7.  0 

0. 500000 

7.  5 

0. 500000 

8.  0 

0. 500000 

8.  025 

0.499781 

8.05 

0.499125 

8.075 

0.498030 

8.  1 

0.496496 

8.  125 

0.494521 

8.  15 

0.492104 

8.  175 

0.489241 

8.  2 

0.485931 

8.225 

0.482171 

8.  25 

0.477956 

(a) 

First  50  cards 

Figure  23.-  Input  data  for  in viscid  calculation. 
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a.  275 

0.473282 

8.  3 

0.468146 

8.  325 

0.462542 

8.  35 

0. 456463 

8.  375 

0.449905 

8.4 

0.442859 

8.425 

0.435319 

8.45 

0.427277 

8.475 

0.418722 

8.  5 

0.409646 

8.  525 

0.400037 

8.  55 

0. 389885 

8.  575 

0. 379176 

8.6 

0. 367897 

8.625 

0. 356032 

8.65 

0. 343565 

8.675 

0. 330479 

8.7 

0.316753 

8.725 

0.302368 

8.  75 

0.287298 

8.  775 

0. 271518 

8.  8 

0. 255000 

8.  81 

0. 255000 

8.  82 

0. 255000 

8.  83 

0. 255000 

8.  84 

0. 255000 

8.  85 

0. 255000 

8.  86 

0. 255000 

8.  87 

0. 255000 

8.  88 

0. 255000 

8.  89 

0.255000 

8.  9 

0. 255000 

8.  92 

0. 255000 

8.  94 

0. 255000 

8.96 

0. 255000 

8.98 

0. 255000 

9.0 

0. 255000 

9.025 

0. 255000 

9.05 

0. 255000 

9.  075 

0. 255000 

9.  1 

0. 255000 

9.  15 

0.255000 

9.  2 

0.  255000 

9.  3 

0.  255000 

9.4 

0.  255000 

9.  5 

0.  255000 

9.  7 

0.  255000 

9.9 

0.  255000 

10.  1 

0.  255000 

10.  3 

0.  255000 

(b) 

Next  50  cards 

Figure  23.-  Continued. 
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10.5 

0. 255000 

10,8 

0. 255000 

11.2 

0. 255000 

11.6 

0. 255000 

12.0 

0. 255000 

61 

31  0 0 

0 

0.0 

0.0 

0.  0 

F 


(c)  Remaining  7 cards. 
Figure  23.-  Concluded. 
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1RMAL  CUOflD* 

STRLTCM  tUP  ALP 

■ 1.300  

J 

AN 

G 

CM 

1 

0.2S10E 

00 

0.0 

0 .77256-04 

2 

0.161  7E 

03 

0 .15456-03 

0.46071-03 

3 

0*63396 

02 

0.766 86-03 

0.13656-02 

4 

0.36 oat 

02 

0 . 19O4E-02 

0.2900E— 02 

S 

0.239  IE 

02 

0.30366-02 

0.5  i486— 02 

6 

0.J720E 

02 

0.64606-02 
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Figure  24.-  Output  for  lnviHcid  calculation. 
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Figure  24.-  Continued. 


120 


AE  DC-T  R-79-4 


ITERATION  HU.  1 


IT 

UP*  AX 

ID 

Jl) 

H HA  X 

JR 

isua 

1 SUP 

RA  VG 

ftp  1 

OF  3 

NS 

1 

0.4026-01 

45 

31 

0*31 2E 

02 

1 

31 

0 

0* 1 50C  QO 

1 *4  00 

0.100 

0 

2 

0*  24  JE-0  i 

4 4 

31 

0 « I62E 

02 

4b 

0 

0.1I9E  00 

1 .400 

0.100 

1 

3 

0.202E-01 

44 

31 

0.1 15E 

02 

45 

31 

0 

Q.I04E  OO 

1*400 

o*  i oo 

1 

4 

0*  1 736-01 

44 

31 

0*1 12E 

32 

46 

31 

0 

Q*  977E-01 

1 *400 

0*100 

5 

5 

0*  1 531;— 0 1 

44 

31 

0.1  02E 

02 

46 

31 

0 

0. 941E-0L 

1 *400 

0.100 

10 

6 

0.  14SE-0  1 

43 

31 

0.95BE 

01 

46 

31 

0 

0.915E-01 

1 »40D 

0*100 

14 

7 

0* 1 J4E-01 

43 

31 

O.904E 

01 

4b 

31 

0 

0*  624E-01 

1.400 

0*  100 

19 

<1 

0. 1236-01 

42 

31 

0.B56E 

Cl 

46 

Jl 

0 

0. 79JE-01 

1 .400 

0*100 

la 

9 

0*  1 17F-OI 

42 

31 

0.CI16E 

01 

46 

31 

0 

0*  7S6E-0 | 

1 . 4 00 

0.100 

21 

10 

0* 1 0 JE-O  1 

42 

31 

0 « 7 79E 

01 

46 

31 

0 

0. 766C-0I 

1 .400 

0*100 

25 

11 

0.S70E-02 

42 

Jl 

0.749E 

01 

46 

31 

0 

0- 752E-01 

I .400 

0*  too 

2a 

12 

0.796E-02 

45 

31 

0 .7236 

01 

46 

31 

0 

0*  697E— 0 1 

1 *400 

0*100 

33 

13 

0. 7d 4c-02 

45 

31 

J*f  03C 

Jl 

46 

Jl 

0 

0-676E-01 

1 .400 

0*100 

37 

14 

0- 7956-02 

42 

29 

0*6a«E 

01 

46 

31 

0 

0*  667E— 01 

1 .400 

0*  100 

39 

lb 

0.7976-02 

42 

3] 

0.o76E 

01 

4b 

31 

0 

0.690E-01 

L .400 

0.100 

47 

lo 

0.  762E-02 

43 

31 

0.660E 

01 

4b 

31 

0 

O.b79£-Ol 

1.400 

0.100 

51 

1 7 

0.&95F-02 

44 

Jl 

0.638E 

01 

46 

31 

0 

0.637E-01 

1 .400 

0*  100 

57 

lb 

O.t-to  6-02 

4© 

Jl 

0*6laE 

01 

46 

31 

0 

0.61 05 '0 1 

1*400 

0*100 

bl 

1* 

0*o5  IE-02 

45 

31 

0*o04E 

01 

4b 

31 

0 

0.599E-01 

1 .400 

0*  LOO 

70 

20 

0*6396-02 

45 

31 

0 .69 |E 

01 

46 

31 

0 

0*  60  IE— 01 

1*400 

0*100 

75 

HtUxa  0 .59E-  01  . CUVH*  0.336-01 


(c)  iteration  history. 
Figure  24. - Continued. 
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1. 125  -0.4362 

43 

0.62  l 

0.358 

0.SV9 

0.899  0.0022 

4 

44 

0.  006 

0.325 

0.  758 

0.754  0.2741 

* 

45 

8.  743 

0.202 

0.561 

0.550  0.6460 

* 

• 

46 

8.000 

0.255 

0.43  d 

0.505  0.851 5 

♦ 

47 

a.  664 

0.25b 

0.450 

0.528  0.8326 

♦ 

4 b 

0.936 

0.255 

0.559 

0.650  0.6404 

4 

49 

9.024 

0.255 

0.63  9 

0.735  0.5035 

♦ 

50 

9.|2o 

0.255 

0.  702 

0.802  0.3837 

ft 

6 1 

9.240 

0.255 

0.754 

0.855  0.2831 

4 

52 

V.397 

0.255 

0.797 

0.899  0.1996 

ft 

53 

9.  504 

0.255 

C.  03  1 

0.933  0.1323 

ft 

54 

9.024 

0.255 

C.  050 

0.959  0.0  809 

♦ 

55 

10. 144 

0.255 

0.877 

0.678  0.0446 

4 

50 

10.592 

0.255 

C.  089 

0-989  0.0216 

4 

57 

1 1.264 

0.25  5 

O.04S 

0.996  0.0088 

4 

50 

1 2.304 

0.255 

0.00  0 

0.999  0.0028 

4 

59 

1 4.624 

0.255 

C.  900 

1.000  O.GOOG 

4 

CO 

2 1.343 

0.255 

0.900 

1.000  0.0000 

ft 

0 ■ 255 

J*  900 

1.000  0.0 

ft 

TOTAL  BODY  JRAj  C0EFF IC  I GNTc  0.21992 


AFTCRBCUY  DMAy  CULFF  tC  ItNf=a  0.16222 


<d) 


Tabulated  solution  and^plot  of 


Figure  24.-  Concluded. 
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TABLE  I.  RELATION  BETWEEN  EXTERNAL  DATA  SETS 
AND  INPUT/OUTPUT  LOGICAL  UNIT  NUMBERS 


Data  File 

{fig  * 16) 

Logical  Unit 

Input 

Output 

1 

14 

9 

2 

13 

8 

3 

12 

2 

4 

15 

10 

5 

11 

3 

123 


AEDC-TR-79-4 


TABLE  II.  BASIC  INPUT  DATA 


ITEM 

NO. 

VARIABLES 

FORMAT 

1 

TITLE  (Three  cards) 

20A4 

2 

NRSTRT , NPRINT , N 3 , ILIM , IITER 

515 

If  NRSTRT  ^ 0 and  N3  f 0,  skip  to  item  12 

If  NRSTRT  f 0 and  N3  ^ 0 stop.  No  further 
input  is  necessary 

3 

LPROG ,N1 ,N2 , IBL , IUNIT,MIT 

615 

4 

GAM, AMINF 

If  LPROG  = 1 skip  to  item  13 

5 

IXY 

15 

If  IXY  = 0,  skip  next  card 

15 

6 

XO,YO  (IXY  cards) 

2F10.0 

7 

IMAX , JMAX , MHALF , KLOSE , LREADP 

515 

8 

DNDZO , XI XM , XM , DS DXIM , XBT , DMAX , XZNEW 

7F10.0 

9 

PLUMIN 

L5 

If  PLUMIN  = FALSE,  skip  next  card 

10 

GAMAP , PTPPFS , AMP , THETAP , TTP , GMP , DELCZ 

7F10.0 

If  LPROG  = -1,  skip  remaining  cards 

11 

LSEP 

L5 

If  LSEP  = FALSE,  skip  next  card 

12 

XSEP , THETS 

2F10.0 

V 

If  NRSTRT  1-  0 stop.  No  further  input  is 
necessary 
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TABLE  II.  (Concluded) 


ITEM 

NO. 

VARIABLES 

FORMAT 

13 

IOPT , K , LVAR1 , LSHAPE , LIC , LDSTAR, LSHPBL 

715 

If  LPROG  / 1 or  if  LPROG  = 1 and  both 

IOPT  = 1 and  LVAR1  = 2,  skip  items  14 
and  15 

14 

NVAR 

15 

15 

XVAR, VAR 

2F10.0 

16 

EL , PT , TT , TWONTT , VISC , RGAS , SCON , DFACT 

8F10.0 

17 

XZ , RLEN , XT , DXP , HLIM 

5F10.0 

If  LSHAPE  0,  skip  next  card 

If  LIC  = 1,  input  CFCl  and  DELTA1 

If  LIC  = 2,  input  CFCl  and  DELSTl 

18 

CFCl , DELTA1  (or  DELSTl) 

2F10. 0 

If  IOPT  = 1,  skip  next  card 

19 

UE1 , DUEDX 

2f 10 . 0 

If  LSHPBL  = 0 and  LPROG  = 0,  skip  items  20 
and  21 

20 

NR 

1615 

21 

XRP , RL 

2F10.0 
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TABLE  III.  SUMMARY  OF  INPUT  DATA  FOR  RESTARTING 


ITEM 

NO. 

VARIABLES 

FORMAT 

1 

TITLE  (Three  cards) 

2QA4 

2 

NRSTRT,NPRINT,N3,ILIM,IITER 

515 

1 

If  N3  = 0 stop.  No  further  cards  are 
necessary 

3 

XSEP ,THETS 

2F10.0 

/ 
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A 

A . 
nD 

A’  . 
nj 

A* 

a 

a 


nD 


B 


n 


n 


C 

C 


D 

f 

G 

g 

H 

I 

J 

Jr 

K 

K 


n 


n 


n 


NOMENCLATURE 

cross  sectional  area  of  inviscid  plume 

coefficients  in  equations  (47),  (48),  (54),  (56),  (60), 
(116)  and  (124) 

coefficients  in  equations  (12,4)  and  (125) 
cross  sectional  area  at  a sonic  throat 
speed  of  sound 

coefficients  in  equations  (61)  to  (63) 

coefficients  in  equations  (47),  (48),  (54),  (56),  (116), 

and  (126) . 

coefficients  in  equations  (61) , (62)  and  (64) 

Chapman- Rube sin  parameter,  equation  (31) 
parameter  defined  by  equation  (82) 
skin-friction  coefficient 
specific  heat  at  constant  pressure 

factor  in  mixing  layer  temperature  profile,  equation  (81) 

weighting  function  defined  by  equation  (39) 

integral  functions  defined  by  equations  (49)  and  (123) 

acceleration  of  gravity 

transformed  shape  factor,  equation  (66) 

integral  functions  defined  by  equations  (50)  and  (58) 

conversion  factor  used  in  equations  (33)  and  (35) 

integral  functions  defined  by  equations  (51) , (52) , and 
(59) 

entrainment  fraction,  equation  £74) 
factor  defined  by  equation  (53) 
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k 

L 

M 


n 


m 

P 

i 

P 

*1 

R 

R 

r 

r 

S 

S 


n 


n 


s 

s. 


T 

T 


U,V 


UB 

Ux 

U,  V 


NOMENCLATURE  (Continued) 

factor  in  equations  (8)  and  (9) 
reference  length 
Mach  number 

1 2 

Mach  number  function  -^(Y-Dm 

functions  defined  by  equations  • (120)  to  (122) 

pressure 

velocity  in  inviscid  exhaust  plume,  equation  (11) 
gas  constant 

functions  defined  by  equations  (117)  to  (119) 
radius 

temperature  function,  equations  (17)  and  (84) 

temperature  function  at  ‘nozzle  exit,  equation  (80) 

temperature  function  in  fully  developed  mixing  layer, 
equation  (83) 

< 

rms  error  of  boundary-layer-edge  velocity,  equation  (5) 

rate  of  change  of  s with  6 , equation  (6) 

rate  of  change  of  s with  x_,  equation  (7) 

, s 

temperature 

constant  in  Sutherland's  viscosity  relation,  equation  (31) 

factor  in  transformed  boundary-layer  equations,  equation 
(26) 

transformed  velocity  components,  equations  (28)  to  (30) 
wake  velocity  factor,  equation  (40) 
friction  velocity,  equation  (41) 

velocity  components  of  physical  flow  field,  figure  4 
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U 

n 


U,  V 


X,Y 


XP 

xs 

■V  JV 

x,y 


a 


y 

<5. 


6. 

i 


Un 


NOMENCLATURE  (Continued) 

velocity  function  at  nozzle  exit,  equations  (75)  and  (76) 

velocity  function  in  fully  developed  mixing  layer, 
equations  (77)  and  (78) 

transformed  velocity  components,  equations  (22)  and  (23) 

transformed  coordinates,  equations  (28)  to  (30) 

coordinates  of  physical  flow  field,  figure  4 

axial  location  of  minimum  ue  downstream  of  boattail 

axial  location  of  separation  point 

transformed  coordinates,  equations  (18)  afid  (19) 

coordinate  of  logarithmic  part  of  boundary  layer, 
equation  (42) 

damping  factor  for  viscid-inviscid  iterations,  equation 
(4) 

eddy-viscosity  factor,  equations  (15) , (16) , (43) , and 
(45) 

value  of  8 for  fully  turbulent  flow 

value  of  6 in  transition  zone  between  laminar  and  turbu- 
lent flow 

ratio  of  specific  heats 
boundary  layer  thickness 

thickness  of  nozzle  internal  boundary  layer  or  thickness 
of  inner  part  of  mixing  layer 

transformed  boundary  layer  thickness 

value  of  5 for  external  boundary  layer  at  nozzle  exit 

f6  u.  r 

boundary  layer  displacement  thickness,  (1  - p|*  )—  <iy 

Jo  Me  e w 
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NOMENCLATURE  (Concluded) 

r6 

■Q 

6*  boundary  layer  displacement  thickness,  (1  - — ) dy 

K J o ue 

e eddy  viscosity  for  mixing  layer,  equation  (86) 

e constant  in  eddy  viscosity  of  mixing  layer,  equation  (86) 

6 angle  ofv  separated  displacement  surface  with  the  axis, 

equation  (2) 

k factor  in  transition. region  eddy  viscosity,  equation  (46) 

y molecular  viscosity 

v y/p 

p density 

t shear  stress 

SUBSCRIPTS 

c refers  to  inviscid  core,  inner  edge  of  mixing  layer 

e refers  to  boundary  layer  or  mixing  layer  external  edge 

I refers  to  inviscid  flow 

0 refers  to  free  stream  reference  conditions 

t refers  to  stagnation  conditions 

V refers  to  viscous  flow 

w refers  to  the  solid  surface  or  the  inviscid  plume  boundary 

x,y  denotes  differentiation  with  respect  to  x,  or  y 
SUPERSCRIPT 

1 denotes  differentiation  with  respect  to  y/6 
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